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I.    IMTPODUCTION 


A.    xi 1 5 TORY 


Hi^h  fr'=qi:ency  ^-a-iio  direction  firdin^,  abbreviated  in 
this  report  as  HFDF,  has  "been  a  topi^  of  interest  since  the 
first  uses  of  radio.  The  science  of  i?DF  has  developed  in 
spurts  of  ingenuity  and  need.  Tne  advances  in  electronic 
devices  a^n.'l  the  more  accurate  jiodeling  of  H?  propa^^ation  have 
teen  important  steps  to  developing  accurate  HIT)?  systems. 
However,  it  has  been  military  necessity  that  nas  spurred  the 
most  important  developments  in  this  field.  The  greatest 
concentration  of  published  literature  on  tnis  subject  resiles 
in  the  technical  reT)ort5  published  durin.^  and  in  tne  decaae 
after  .^orld  War  Two. 

The  majo^  studies  of  H7D?  have  dealt  with  the  problems  cf 
polarization  of  skywave  signals,  the  effects  of  multipatn 
propag'=' tion ,  the  statistics  of  H?  propagation  through  the 
icncsphere  and  the  ^^  e veloDmer  t  of  rlFDF  antennas  and  arrays  of 
antennas.  A.n  iTiportant  distinction  nas  developed  from  these 
studies.  There  are  three  tyoes  of  HFD?  antennas:  (1,.  wiae 
aperture.  '' 2  ■  medium  aperture  and  (3,^  narrow  aperture 
antennas.  .vhether  an  antenna  array  cc^oares  pnase  or 
amplitude,  the  p^imar^y  neasu.re  rf  problems  attendant  to  its 
accurate  i^pe-^aticn  is  the  width  of   its   aperture.   Tne   term 
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aperture  in  this  paper  refers  to  the  linear  spatial  extent  of 
an  antenna,  not  to  an  area.  The  unit  of  r.easurernent  is  either 
meter  or  wavelength.  If  the  aperture  is  on  the  order  of  O'-^e 
qua'^ter  c :*  less  of  a  wavelenrtn.  it  can  te  oonsider^ed  narT'ow 
aperture.  and  it  suffers  the  greatest  rumter  of  difficulties 
to  achievin,^  accurate  direction  findin>-;  capability. 

During  '.vorld  War  Two  the  Allies  experienced  considerable 
success  v/ith  the  landbased  medium  and  wide  aperture  systems 
and,  not  unejpe ctedly .  limited  success  with  narrow  aperture, 
shipboard  'd7d?  systems.  One  of  the  major  tar.'ets  of  the 
shipboard  H?j?  systems  was  German  submarines.  As  improvements 


were  made  to  all   tvoes 


HFD?   systems,   tne   submarine's 


transmissions  became  more  and  more  vulnerable.  In  an  effort 
to  maintain  communications  and  to  tnv/art  'd?D?  systems,  t.ie 
Germans  shortened  the  duration  of  transmissions  to  lower  tne 
probability  that  the  transmissions  would  be  intercepteo  and 
SLbsequentlv  located  by  "JD?  .  A  highly  effective  means  of 
shortening  transmission  time  was  to  record  the  information  on 
tape  and  then  to  transmit  via  the  radio  at  a  much  faster 
playbaC'^  SDeed  .  When  this  metnod  was  coupled  with  the 
practice  of  economizing  on  the  amount  of  inforKiation  sent, 
signal  durations  were  shortened  by  more  than  an  order  of 
magnitude.  A  U-bcat  employing  such  measures  was  appreciably 
less  susceptible  to  HFD?. 

The   problem   of  locating  a  short  duration  signal  remains 
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today.  It  is  still  ^  ccmrion  p^otlem.  even  vhen  t'ae  ta^r^et 
t^ansritter  dees  '^ot  attempt  to  corrpT'ess  its  sir:;iial.  In  a 
tactical  situation  it  is  typical  that  tne  conrrunica  ticn  net 
control  station,  usually  co-located  with  t.ie 
of  f  icer-in-tac t ical-conmanc,  will  act  as  a  broadcas  ti^ig 
statioTi.  and  the  outstations  will  not  transmit  or  will  only 
transmit  a  t'^^ief  sipnal.  In  the  case  cf  'nanuai  Tiorse  this 
signal  ray  be  an  "^"  for  "rcf^er  -rour  last  transnissicr",  cr 
in  the  case  of  tactical  voice  c  or^nunica  t  i  c^s  tne  outstati^n 
will  trieflv  key  the  ricrophore.  In  eitner  case  tne 
transmission  r\ay  not  last  lonrTe^'  than  200  to  4<?2 
milliseconds. 

The  rapid  growth  of  diijital  c  onrun  i cat  ions  nas 
significantly  increased  the  ease  with  whicn  a  curst  signal 
can  be  generated  and  " r el iabl^'-'"  received,  ^^-iven  a  digital 
pi'lse  of  duration  "t"  and  a  total  signal  duration  of  "t", 
there  is  a  si'^Dle  ej^ression  for  the  amount  of  information  in 
bits  that  can  oe  transmitted. 
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figure  1 
Partial  short  duration  signal  bit  stream 
If   the   number    of   bits   is   "n",   tnen   N=T/2t. 


IS 
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approTimatelv  the  reciprocal  of  the  ■bandwidth;  tnerefore, 
N  =  (T)(3'*).  a'Dprcxi'^atelY.  In  tne  high  freouencv  range  a 
bandwidth  cf  12  .-cHz  can  "be  readily  achieved.  As  an  example. 
if  BW=ir  khz  and  T=5?Z  pis.  N  equals  52''?2'  bits.  Five  tacusa^d 
bits  is  sufficient  to  provide  ccrsiderable  encoded  tactical 
irf  crna  t  icn .  "ven  if  tne  si^^nal  diiraticn  were  linited  to  2j0 
[PS  and  the  ^i^st  half  of  the  signal  dedicated  tc  alerting  tne 
destination  receive^,  the'^e  vonld  ren^ain  10Z0  bits  for 
inf  orr.a  t  ion .  At  five  bits  pe^  si'^tol  and  an  average  of  five 
symbols  ve^  word.  this  would  allow  forty  words  to  be 
communicated  in  the  sDace  of  2?d  ms.  ^^'or  a  shin  on  a  covert 
thirty  day  patrol  sending  a  single  daily  status  report  via 
burst  corrmr  nicat  ions  .  the  total  communications  trans~iis  si  on 
tine  would  amount  t'^  six  seconds  (about  two  one  milliont.is  cf 
tne  pa  t  '"n  1  perio^i  )  . 

To  date  this  author  has  not  been  able  to  find  any  past  cr 
current  research  on  the  capacilities  of  H?L?  systems  to 
exploit  short  duration  signals.  This  is  probably  due  to 
several  reasons.  The  primary  reason  is  tnat  H?D?  engineers 
are  still  absorbed  in  the  more  basic  problem  of  improving 
HFDF  against  medium  and  long  duration  skywave  signals. 
"S^speciellv  in  the  case  of  narrow  aperture  hY'D?  antenna 
systems,  there  remains  considerable  need  for  improvement.  In 
the  case  of  wide  aperture  '-^?D?  antenna  systems,  tne  problem 
of  short  duration  signals  seems   to   be   tractable.   however. 
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there  does  not  appear  to  exist  a  comprenens i ve  stady 
addressing  this  proclen,  and  to  conduct  such  a  controlled 
study  of  skyvave  propagated  signals  would  ce  expensive. 
Notvit hs 1 5nding  these  dif f icrl t ies  .  the  short  duration  signal 
CGuld  easily  oeccne  an  acute  tactical  DroV.lem  for  the  side 
that  carnot  exT)lcit  it,  and  the  prohleT,  tnerefore  deserves 
i^imediate  attention. 

B.  FURPCSZ  CF  THESIS 

The  general  question  '^f  inte'^est  is  how  good  are  existing 
h?DF  syster-s  at  determining  lines  of  bearing  on  sno^^t 
duration  signals.  Any  conplete  examination  of  an  dFD?  system 
requires  one  to  investigate  tne  characteristics  of  either 
wide  0^  narrow  aperture  antennas.  Additionally,  one  r.ust 
exanine  the  p^^ohlens  of  site  location.  signal  acquisition 
syster,  "eceiver  d-^nodul  ^  tion  ,  tei'^in^  sense  circuits  and  tne 
noise  en  vi  ^'onner  t  .  Perfnrnance  ni-st  also  "be  deternined  f^^ 
ground  versus  skywave  and  rnultipath  versus  single  patn.  rnis 
would  "be  an  enormous  task:  if  all  current  HFDF  systems  we'^e 
considered.  The  scope  of  this  effort  is  much  mo^^e  restricted. 
The  purpose  of  this  thesis  is  to  investigate  the  mean  and 
variance  of  tearin.-?  estimates  for  short  duration  s'-cywave 
signals  received  with  a  narrow  aperture  HFDF  system.  Several 
statistical   procedure's   are   developed  and  compared  with  the 
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intention  of  li  scrininati^^  between  reliaMe  and  unreliable 
data  and  calculating  the  best  bearing  estimate  fron  tne 
reliable  data. 

The  problen  of  hov  to  acquire  a  snort  duration  si-^nal  or 
how  to  interface  such  an  acquisitio'^  syste-n  to  an  HPDi'  system 
will  not  be  addressed  in  this  report. 
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II.  TH"OPiTICAL  COiSIDZPATIONS 


A.  hi;. I  JPEQUZMCY  SKYV/AV5  CHAMNEL 


Hi^h  frequency  '?.?^  is  the  re,^ion  rron  three  to  tnirj 
regahertz.  Its  corpleTitv  is  due  to  '^anv  latural  pnencneria 
which  are  in  terdep^n*^  ent  ,  co;~plex  in  thensslves  .  so"ie  pooriT 
understood,  cosric  and  ricT^ocosni?  in  extent  and  diffici;lt  to 
measu'^e.  The  p^ina^v  oonpleritv  fo"  the  '"I?D?  en^^ineer  resides 
in  the  ionosphere.  The  ionosphere  can  be  considered  an 
inhomOf-eneous  plas-na  trxat  surrounds  a  sphere  of  finite,  jut 
variable,  cond^ict  iv  i  tj  and  separates  the  sphere  fror,  free 
space  and  nMT«ero!TS  solar  and  galactic  sources  of 
disturbances.  The  ionosphere  has  been,  and  continues  to  be,  a 
subject  of  considerable  research.  References  1  tarou^^n  '-  are 
a  rich  "'^esoi'i'''ce  of  in^'^orr'ation  en  the  io'^osr^nere  and  "esearch 
repoT-ts  are  added  monthly,  but  the  scnpe  of  the  problem  is 
imrense.  To  predict  accuratelv  ionospheric  corcitions,  one 
r^ust  be  equipped  with  ro^e  than  the  phvsics  of  the 
ionosphere.  The  physics  provides  the  equations  of  tne  syster, 
but  the  forcinr:  functions  and  the  boundary  conditions  nust  be 
sufficiently  measured  to  forecast  accurately. 

The  fcrcin;-^  functions  are  tne  solar  flux,  gra  v  i  t  a  li  onal 
waves,  weather  related  avnapics  and  the  si.-^nal  of  interest. 
The  ionosphere   is   usually  ncdeled   in   terms   of   electron 
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concen  t  ra  ti ''^p.s  *  the^efce.  phenonera  v/hich  affect  tne 
concentration  o^*  the  excitation  ^f  the  plectrons  irive  tie 
system.  The  solar  flux  is  primarily  a  diurnal  paenonenon;  its 
impact  is  strongest  in  tne  portion  of  the  plasma  ill-jminated 
by  the  si:n.  This  fl^.^x  is  naae  up  of  elec  troma.-^net  ic  energy 
and  streams  of  particles.  In  the  case  of  sun  ST^cts  and  solar 
flares,  there  are  often  increased  emissions  tnat  tend  to 
disturb  the  normal  structure  of  the  ionosphere,  f Perhaps  it 
would  he  mc^e  accurate  to  state  that  tne  disturoance  is  to 
our  model  of  the  ionosphere.'  Tne  solar  disturbance  evcl^'^es 
in  three  stages.  The  first  is  the  impact  of  elect  ^oma.^net  ic 
energy  in  tne  ult-^aviolet  and  x-^^ay  '^a'^.^^es  t  na  t  causes  an 
increased  electron  concent-cation  in  the  lovest  electron  layer 
(D-layer^.  The  second  effect  is  the  arrival  of  hi?;h  ener^.-y 
protons  and  alpha  particles  that  also  increase  tne  D-laye". 
The  duration  of  the  disruptions  due  to  these  tv/o  "oaases  is 
limited  to  several  hours.  1!r^e  third  Dhase  is  lie  arrival  of 
low  ener,2-y  protons  and  electrons  whicn  shower  the  earth  in 
patterns  molded  b"""  the  Earth's  magnetic  field.  In  this  pnase. 
which  may  last  as  Ic^r.  as  several  days,  the  ionosuhere 
experiences  magnetic  storms,  an  increase  in  the  D-layer. 
sporadic  conditions  in  the  next  higher  ^'-layer  and  tne 
spectacular  au'^o^'ae. 

Acous  to-^ravi  ty  waves  constitute  a  forcin,^^  lunction  of   a 
different   sco^e.   Periodic  variations  in  tne  i/namics  of  the 
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Earth-nccn-sun  gravitational  systen  and  isolated.  anoinalo'^s 
gravitational  activity  en  the  Earth  con"bin9d  with  Hf  acoustic 
waves  (The  Mt .  St.  Helen  eru"ption  vas  a  recent  source  of  S'lch 
waves)  exert  forces  that  distort  the  .^ene^^al  concent'*!  j 
spherical  form  of  the  electron  plasma  layers.  The  distortions 
are  not  only  static.  Travelin..:?  ionospheric  disturbances  a^e 
not  uncorrron.  and  their  effect  is  to  create  a  doppler  shift 
on  transmitted  signals.  If  the  ionospheric  disturbance  is 
tilted,  the  ray  trace  of  a  transmitted  si^:nal  will  ce  bent  in 
azimuth. 

The  third  mentioned  lorcin.^  function  is  the  weatner.  '^■i.e 
dynamics  of  the  weathe^  affect  the  pressure,  the  temperature 
and  the  mlTint:  of  the  atmosDhere.  These  th'^ee  factors  in  tu?"n 
have  a  sii:;ni  f  i  cant  impact  on  the  elect^-on  ccncen  tr-a  ti  ons  , 
particularly  the  concentrations  at  the  lower  altitudes.  Tne 
weather  is  also  a  veT'y  important  factor  in  ni/;n  frequency 
ranges  because  it  is  a  noise  source.  '"luch  of  tne  hiijh 
frequency  bacfc^round  noise  is  attributed  tc  thunderstorm 
activity  which  is  cent  innall;?'  occurrin-?;  at  some  point  on  tne 
Earth.  (It  should  b^  noted  that  ^^nst  of  the  elect  rcma^net  ic 
enerf:y  cf  a  thunder's  term  is  in  the  VL?  ■^^et-:;!  on , ) 

^anmade  signals  a^e  one  of  the  smallest  forcings:  functiors 
acting  on  tne  ionosphere,  but  they  are  naturally  of  ^reat 
interest.  The  tarsret  si^-nal  iniects  itself  into  the 
ionosphere;  it  operates  on  the  ionosphere  and  is  operated  on 
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ty  the  ionosphere.  The  st^idy  of  this  interaction  has  lead  to 
a  description  of  the  change  of  transnitter  anten-.a 
polari?ation  to  elliptical  polarization,  the  phenorena  of 
refracted  hi^h  freonency  v/aves,  nultipatn  interference  and 
the  conceTits  of  maTiTium  useable  frequence'-  (MUF),  lowest 
useable  freauency  (LH?)  and  optinum  working;  frequency  (TOP'. 
The   ionosphere   is   a   systen   vitn   largely   fluctuating^ 


boundary  conditions.   The   sur 


of  the  Zarth  is  tiie  only 


boundary  that  can  be  considered  fixed  with  respect  to  daily, 
seasonal  and  eleven  year  solar  cycles.  Other  boundary 
conditions  are  much  ncre  dynamic.  Of  tnese.  tne  layering  of 
electron  concentrations  is  priTiarv.  The  inner  two  layers  v  J 
and  I,  which  are  -ncstlv  the  result  of  solar  electronagne  t  i  c 
radiation  have  been  mentioned.  Tne  cuter  layer.  7,  v/aicn 
often  is  subdivided  into  an  ?1  and  "2  layer  is  relatively 
r  0  r  e  stable.  It  ^  e  "ia  i  n  s  when  the  "o  c  r  t  i  o  n  of  the  i  o  n  c  s  p  n  e  r  e  of 
interest  rotates  into  the  sola"  unbra  and  tne  D  and  !■"  lavers 
disperse.  The  0  and  E  la-'-e^s  du^'ing  davlignt  a^e  responsible 
for  the  non-deviat i ve  attenuation  of  "uch  of  tne  nF  spectrum 
of  interest  (3-12  MHz  \  The  dispersion  of  the  D  and  E  layers 
permits  the  ?  layer  to  become  a  virtual  reflector  situated  at 
altitudes  typically  from  220  to  ^-dd  Km.  iTne  actual  mechanism 
of  ■orcpa;?at  i  on  through  the  "^  layer  is  refrection  whicn  can  be 
modeled  as  reflection  from  a  virtual  heignt  greater  t  lan  the 
actual  zenith  of  the  bending  ray.'>  Y   T5ror)agat  ion  opens  up  tne 
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dv.e   solely   to    the 


evening;  airwavs  to  lor..^  distance  cc^mur  ica  t  ions  and  attendaiit 
long  distance  ?.?D?  in  the  3-12  MHz  '^ane^e .  7or  tne  en<?ineer 
this  is  a  rixed  blessing. 

The  ri?D^  engineer's  interest  in  s  cywave  propagation  is  in 
the  difference  "between  the  direction  of  arrival  of  the  tdr,-:^t 
sij^nal  and  the  g^ea^^  circle  bear  in.-;  to  the  tar,^et  and  in  tne 
variance  of  the  reas^jrernen t s  of  the  ar^Tle  of  arrival.  Aside 
f^'OiT  the  eq^jip^ert  liritaticns  and  site  lo:ation  distortior.s 
and  "ef  lect  i  0'"  5 ,  r^ny  of  the  errors  and  variances  that  need 
to  he  "esolved  to  irp^^ve 
ionosphere. 

In  the  evening,  targets  of  interest  in  the  3-12  MHz  hand 
can  he  exnlcited,  hut  there  is  a  considerahl/  greater  chance 
of  in t er f eT'ence  fron  other  discrete  sources  or  fror^.  general 
noise  sources.  Additionally,  tnere  is  increased  complexity 
when  signals  routinely  arrive  after  two  cr  three  hops  v.'hicn 
correspond  to  naxinun  distances  of  £003  and  12.^'^?  --^i . 
respectively.  Over  these  distances  the  errors  and  variances 
due  to  in  te^ef  eT'ence ,  fading,  tiltirg  a'^d  scatterino"  increase 
to  a  point  that  even  vide  aperture  antennas  cannot  produce 
useful  fi.T  information. 

An  inpcrtant  consideration  for  narrow  a^^erture  antennas  is 
that  the  effective  anerture  of  one-qusrter  v/avelengtn  at  <:d 
MHz,  a  typical  longhaul  davtiie  frequency,  oecones  a 
cne-sixteentn  wavelength  at  5  MHz.  a  nigntti-ne  frequency.  I'ne 
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loss  of  effective  ^"nerture  further  exaceroates  tne  prcblem  cf 
dete^ning  a  bearing  ar.d  its  variance.  The  effective  hei-^ht  of 
the  arter.na  is  also  a  function  of  frequency?  tnerefore,  ore 
can  expect  the  array  patte^^n  to  chan^^^e  vith  txie  change  in 
operating  frequencies. 

B.  NARRCW  A?i?.TUR^  D?  AMr^NNAS 

The  knowledge  of  the  ionospnere  has  grown  extensively  in 
the  past  forty  years.  Investigators  can  nov/  feel  reascnaoly 
confo^tacle  vith  the  developed  models  and  the  in^proved 
sensors.  especially  the  ex  t  "a-teT-^'es  t  rially  sensing 
satellites.  leneral  predictions  are  possible  and  a  new 
favorite   c^nnpiiter   aid   is   the    software    that    predicts 


propagation   and   disDlays   rav   t'^acin^^s   (see 


r  i:  1    . 
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Appendix  C'.  An  'I?"!??  engineer  can  review  the  general 
propagation  scheme  with  an.  assurance  that  he  unierstands 
sufficiently  tne  problems  presented  by  the  ionosphere.  Sut  in 
the  case  of  narrow  aperture  r.FD?  antennas  one  niust  guard 
against  the  feeling  of  confidence  induced  by  a  knowledge  cf 
the  general  situation.  One  is  reminded  of  the  situation  where 
a  blind  man  feelin.-"  the  trunk  of  an  elepnant  attempts  a 
5;eneral  description  of  the  elephant.  In  the  case  of  a 
1.5-meter  aperture  antenna  sampling  a  wavef^ont  in  tne 
60-meter   band,  the  dimensional  com"oarison  with  a  nand  and  an 
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elephant  is  accurate. 

There  are  two  comrrionly  used  ty^es  cf  narrow  aperture 
antennas.  One  type  relies  on  amplitude  conDarison  to 
determine  direction  of  arrival  and  the  second  type  relies  en 
phase  corpa^ison.  An  exap^ple  of  tne  fcr-^er  is  tne  si-nple  lc,-p 
and  of  the  latter  is  the  Adcocn:  .  (Reference  5  points  out  tnet 
the  phase  and  amplitude  dis ti n? ti c  is  not  clear  cut  in  the 
case   of   the   Adcock.''   The   case   of   the   sir^ole   loco   is 


illustrated  Delcw. 
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Figure  2 
SiFTDle  loop  sensing  direction  of  arrival 

The  direction  of  arrival  of  the  signal  is  determined  hy  ttie 
T'elative  orientation  ^f  the  loop  and  the  horizontal  component 
of  the  ragnetic  field. 
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The  case  of  the  Adcock:  (actuallv  one  half  of  a  U-Aijoc'r^)   is 
illustrated  as: 


^ 


^ 


V 


r~-  -  -^~, 


"-^^ 


8  0^   s\h'^  (vk") 


'Fii??ure  3 
Simple  Adcock  sensing  direction  of  arrival 


The  diT'pction   of  arrival  of  the  sisr.al  is  deternined  ty  t.ie 
phase  difference  "between  the  ^vo  elenerts. 

The  two  examples  above  only  serve  to  illustrate  now 
direction  of  arrival  information  is  deterrrined.  Real  syste-^s 
employ  more  elements  to  resolve  amhi.'^ui  ties ,  improve  acciiracy 
and  enhance  resistance  to  noise.  The  point  is  tnat  trie 
fundamental  process  relies  on  an  element  sensing  amplitude  or 
phase.  This  fundamental  rrooess  is  in  tu^n  the  fundamental 
difficulty  f'or  narrow  aperture  HFD^  antennas. 

In  Kef.  6  Gethin.^  uses  computer  simulation  to  plot  w^ve 
interference  of  multimode  si^^nals  in  terms  of  surfaces  of 
constant  phase  (CP)  and  constant  amplitude  (CA),  Sot  the 
ideal  case  of  a  single  specular  component  v.ith  no  scat  teri  n.-; , 
the  surface  cf  C?  and  CA  is  a  plane  whose  normal  is  the 
direction  of  oropagation.   In   the   case   of   two   rays.   t;ie 
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interference  patte"^n5  '^ep^esentei  oy  the  surfaces  of  0?  ard 
CA  vary  with  the  anrular  separation  of  the  rays  in  elevation 
and  azimuth.  In  all  of  the  patterns  p'^esented  in  F.ef.  5  in 
which  the  anplitr^des  of  the  two  conponent  rays  differ  by  only 
ten  percent,  "lajor  distortions  to  the  ideal  planes  of  CA  ar^d 
CP  occur.  ApproTirat  el  V  plana''  portions  of  the  surfaces  of  JA 
and  C?  extending  to  se^^eral  wavi=len<?ths  in  len.-rth  are  up  to 
sixty  de.^rees  diffe^^ent  from  the  t'^ue  an^rle  of  arrival.  Tnere 
are  also  -cin'^s  in  the  uhase  fT'orts  that  vary  tne  phase  up  to 
ninety  decrees  in  less  than  the  space  of  one  vavelene:th.  In 
the  cases  where  rce  than  tv/o  ^ays  are  present,  the 
interference  patterns  ■becon'=  ruch  po^e  conplex. 

It  is  obvious  that  the  spatial  extent  of  v/ide  aperture 
antennas  is  needed  to  resolve  sucn  interference  patterns  in  a 
short  perici  of  time.  Balser  and  Smith  in  Ref.  7  explained 
that  when  the  outputs  of  two  antennas  were  correlate!,  the 
antennas  had  to  he  spaced  forty  wavelen^^tns  in  tne  case  of 
sin<gle  hop  and  ten  wavelengths  in  the  case  of  -nultinop  to 
lower  the  correlation  coef f ici'^^nt  to  -^.5  .  For  a  narrow 
aperture  antenna  to  detect  phase  or  amplitude  distortions  of 
this  naenitiide.  the  tine  of  observation  rust  he  relatively 
long.  However,  it  is  necessary  to  detect  such  distortions  to 
permit  an  assessment  of  reliabilty  to  he  assigned  to  bearings 
measured  in  distorted  fields. 

The  description  above  of   interference  patterns  was  for  two 
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sources  nearly  equal  in  amDlitude.  T'ne  condition  of 
comparable  a'^i'li  tudes  is  one  that  res^ilts  in  severe 
distortion.  As  the  amplitude  of  one  cf  tne  rays  becomes 
substantial!?/-  less  than  the  other  ^ay.  the  interference 
pattern  approaches  the  ideal,  undistcrted  planar  pattern 
(implicit  is  tne  assumption  tnat  single  mode  scatter  is  also 
very  weak:).  Assumin,^  that  at  least  one  cf  the  rays  of  a  two 
ray  interference  -oattern  is  fadinr.  it  can  be  expected  ta<3t 
for  short  periods  of  time  ideal  planar  CA  and  C"  waveironts 
can  be  obse^^ved.  There  is  no  readv  means  of  identifyins  these 
moments;  however,  if  the  fading  is  random,  the  planar  CA  and 
C?  wavef rents  should  1°  the  statistical  mean  of  tne  measurea 
wavefronts.  The  ^ate  of  fading  sh'^uld  therefore  be  a 
parameter  to  indicate  the  time  duration  required  to 
statistically  acquire  a  measurenent  of  the  true  anj:;;le  -^f 
arrival . 

It  is  noted  in  Ref.  6  that  for  a  single  ray  with  i'araday 
rotation  induced  elliptical  polarization  the  fade  rates  are 
measured  in  seconds  per  cycle.  In  ?ef.  12  polarization  fading 
with  periods  of  10  seconds  snd  2d  db  fade  deptns  were 
reported  as  common.  If  two  or  more  rays  are  present,  fading 
is  measured  in  cycles  per  second.  This  indicates  that  a 
narrow  aperture  H?D?  system  will  reauire  appr  ox  iima  teiy  a 
second  to  recognize  the  fading  condition  if  strong  multipath 
interference    exists.   The   time   reauired   to   average   the 
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interference  "oattern  is  related  to  the  polari  zat  ioi  fadin^^  oi 
tne  doninate  Tiode.  The  amplitude  of  a  polarization  faded 
signal  is  a  stocnastic  process;  therefore,  tnere  is  no 
deterministic  functional  rela  ti  onshiT)  "cetween  time  and 
fadin:?.  A  measure  of  the  rapidity  vith  wnich  facing  is 
fluctuating  can  he  expressed  in  terms  of  a  f  ad  in.5  power 
spectrum.  (Section  5. 4. 3  of  Ref.  1  discusses  tne  concept  of 
fading  Dover  SDectrr.m.)  If  tnere  is  a  lar^'e  portion  of  tne 
"fadin.^  Dcwer"  in  the  higher  frequencies  {13^!  to  12^6  Hz  • , 
the  faain.^  is  fast.  If  the  "fading  power"  is  primarily  in  tne 
the  0.1  to  1.2  Hz  region,  the  fading  is  slew.  In  tne  case  of 
polarization  fading  it  nas  already  been  noted  tnat  fading  is 
typically  in  the  seconds  pe^^  cycle  range.  Therefore,  an 
antenna  whicn  does  not  have  sufficient  spatial  aperture  to 
average  interference  patte^-ns  must  rely  on  fading  to  permit 
time  averaging.  The  time  required  for  averaging  is  a  function 
of  the  interference  pattern  and  appears  tc  be  on  the  oraer  of 
five  to  ten  seconds. 

A  measurement  experiment  reported  ty  iain  in  Kef.  8 
demonstrated  how  time  averaging  of  hearings  reduced  the 
variance  associated  with  the  mean  heaT'ing.  Using  a  U-AdcoctC 
with  huried  feeders.  bearings  on  skywave  signals  were 
recorded  at  five  bearings  per  second.  An  autocorrelation  -^'f 
the  bearings  was  computed  and  the  resulting  curve  was 
approximated  by  the  exponential  expression; 
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R(  T  )    =   exvi-  T'/To 


wnere  is    a    r)arameter   associated   with   fitting    an 

exponential      curve      to      the      neasured      "bearings.      The    fornula 

^      2. 

relating  the  variance  of   the   ^lean   oearing   (Otj   and   tne 

variance  of  a  single  observation  (C)  is: 
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where  T  is  the  tine  interval  over  which  the  hearings  were 
averaged.  5ain  reported  that  for  {%":  =  0.56  (corresponding  to 
considerable  hearing  fluctuation),  the  variance  was  reduced 
by  a  factor  of  12!  in  12  seconds.  The  12  second  duration 
roughly  cc^^esponds  in  order  of  magnitude  to  the  reciprocal 
of  an  average  fade  T'^te. 

G.  SUM'^'ARY  0?  THEORETICAL  CCN SI  DERATIONS 


The  narrow  aperture  H?D?  antenna  is  physically  linited  to 
tire  averaging  operation  against  s^ywave  signals.  In  tne  case 
of  ideal  ionospheric  propagation,  the  antenna  system  can 
perform  within  equioment  and  site  limitations.  If  tne  site 
errors  are  ■^'^own,  the  equipment  and  array  calibrated  and 
tnere  is  a  good  SMR,  average  bearing  errors  cf  i).o  to  1 
degree  and  variance  of  5  degrees  squared  should  be  possible. 
If  multipatn  propagation  exists  with  fading  up  to  Zi    db.  time 
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aver3s;ing  over  at  least  ten  seconds  with  sampling;  at  a  tout 
five  T)er  second  should  reduce  iiost  of  the  variance  due  to  tae 
complex  interf e^^ence  patterns. 

The  difficulty  of  ottainin^"  accurate  Hi'D?  a^^-ainst  short 
duration  signals  usinf:  a  narrow  aperti:re  array  is 
considerable.  In  the  case  of  n u  1 1 i p a t h  Interference  in  v; h i c h 
at  least  tv/o  rays  are  co-parahle  in  anplitude,  the  l)F  error 
on  a  short  duration  signal  witn  only  one  sample  bearing  could 
te  up  to  ninety  degrees.  This  is  the  extreme  of  tearing  error 
due  to  phase  and  anplitude  front  distortion.  In  a  none 
hospitatle  multipath  environment  the  system  performance 
should  te  much  better,  but  there  is  little  experimental 
evidence  by  which  ere  can  assign  typical  tearing  errors  ard 
variances.  The  analysis  of  the  narrow  aperture  antenna  system 
in  the  following  sections  provides  Derf or'~'ance  data  on  ar. 
experimental,  state  of  the  art  system. 
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III.  SVRI  H~LF  ANT5MNA  SYSTEM 


A.  INTRGCUCTION 

The  ?lec t r oragnet ics  Division  of  Southwest  Pesearch 
Institute  ''S'-'/RI),  located  in  San  Antonio.  Texas.  nas 
developed  and  tested  a  nev;  desi-;n  for  a  narrow  aperture  r.FDF 
antenna  s'^sten  to  operate  a.?ainst  "both  ^^round  va^e  and 
Sivyvave  si-^nals.  The  significance  of  tnis  nev;  aesi^n  is  tnat 
it  is  a  -^ast  nountahle  narrow  aperture  antenna  that  is  a 
fixed  arrav.  There  are  no  novin.^  parts;  tnerefore,  it  is 
ideal  for  the  shipboard  enviror~'ent.  T'-^e  ele-^ents  of  tne 
array  are  si'^ple  loops  and  spaced  loops.  Tne  latter  will  be 
shown  to  have  polarization  independent  qualities  and, 
theref^^e.  to  0"=  ideal  fo^  erploitin?  s'^-'-waves.  The  array  and 
the  associated  i ns trunen ta t ion  of  the  system  are  hi^h  speed 
and  ccr^puter  controllable. 


T  n 


he  prira^y  ^efe'-ence  for  the  analysis  that  follows  i 


.n 


in-house  repoT't  prepared  by  the  syster  architects  [itef.  9  >  . 
The  system  herein  described  has  beer  patented.  The  autnor  of 
this  thesis  has  visited  the  San  Antonio  site  ini  has  operated 
the  H?D7  system  wit^.  the  assistance  of  the  Sv/?.I  "oe^sonnel. 
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3.  THEORY 


To   understand   the  operation  of  the  spaced  loop  arrav  one 


nust  re^riew  the  thecrv  of  the  si^^ple  Iogt).  The  figure  on 


the 


next  pa.<^e  deoicts  a  sinple  Iood  set  i  ^^  a  coordinate  s/sterr 
with  an  incomin.^  si<?nal  ray  (Z  fielc'  cor^ponents  laceied  Zv 
ani  ^h)  .  The  an;?le  pni  (<$))  is  tie  azin^utn.  measured  in  tne  XY 
plane  of  the  in coming  ray.  The  plane  of  the  loop  is  aliened 
with  the  IZ  plane.  The  incidence  an.gle  tneta  {©^  is  measured 
in  the  plane  defined  ly  the  Z  axis  and  the  si^:nal  ray.  The 
signal  is  considered  to  have  a  vertical  and  Horizontal 
electric  field  component  (Ev  and  Ih ) .  The  expression  for  the 
output  "clta^e  (sinusoidal  input,  time  variation  suppressed,' 
is  : 

VI  =  -  ( 1 V  '  ^  CC  5  4>  -^  f  ^h  ' )  (  COS  e  5  i:;  ()j.  )  exD  i  j  4>  h  '•    (1) 

where: 

iv ' :  relative  amplitude  of  the  vertical  component 
Fh ' :  relative  amplitude  of  the  horizontal  component 
4>  :  a  zinuth 

B  :    angle  of  incidence 
yh:  phase  of  horizontal  component  relative  to  tie 
vertical  component 
VI:  simple  loop  output  voltage 


With  respect  to  sk-.-vave  signals,  tne  significance  of  equation 
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(1)  is  that  it  is  polarization  leperdeit.  Hi'Dj?  systens 
(generally  "ely  on  isclatin,?  a  null  in  the  array  pattern  tnat 
can  be  related  to  the  azinuth  of  tne  incnnin^;  si^-^nal.  Tne 
null  used  should  onl;''"  "be  a  function  nf  the  ta''£:et's  bearincT. 
The  sirrple  loop  works  well  with  .^^ound  waves  fo-  waich  case 
theta  is  equal  to  93  degrees.  When  theta  is  90  decrees 
equation  ''1^  reduces  to  a  function  of  only  one  spatial 
variable,  phi.  vhicn  is  the  desired  bearing.  The  output 
voltage  in  this  case  is: 

VI  =  -Sv'  COS  <t> 


The  sinple  loop  does  not  function  acceptably  a^^ainst 
skyvaves.  In  the  case  of  skyv;aves  the  loop  vcltaf:;e  is  a 
function  of  the  two  spatial  variables,  theta  and  phi.  and  t::e 
relative  phase.  The  nulls  coated  bv  these  three  variaoles 
are  too  nurarous  and  the  available  ii^easurer^en  ts  too  sparse  to 
resolve  all  the  arbiguities. 

A  solution  to  tne  polarization  dependence  linitation  of 
the  si '^1  pie  loop  is  to  ccnoi'^. e  two  si'nple  loops  into  a  t'^o 
element  interferometer  as  illustrated  in  fi,^ure  2.  The  loops 
are  connected  in  parallel  with  o"OT}osin,-^  phase.  Tne  output 
voltage  of  t-iis  array,  i-cnovn  as  coaxial  spaced  Icops,  can  be 
deterrined  bv  Dattern  mul t iDlicat ion  [Pef  .  IZ  ].  The  pattern 
of  I h e  spaced  loop  Fiv^^y  is  equal  to  the  p r o di' c t  of  tne  .^ r o u p 
pattern   and   the   elerent   pattern.  The  .^-^oup  pattern  cf  the 


Figure  5 
Coaxial  Spaced  Loops  and  Coordinate  System 
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array    is: 


Gr    =    J    0    d    SIN(  e  )3I\(  Cj)  ) 

d :    separation    cf    the    tvc    loops 


/3:    2fr/7s 

A*    si^'nal   wavelenf=;th 

Nov   let: 

Ev    =    j2v'      d/2 
in    =    2^^2.'      d/2 

This  perr^its  the  spaced  loop  output  voltage  to  le  yritten  as 

Va  =  -Ev  SIMSSIN  2(^-  Eh(SI\'  203  I  M  4>' P7p '' J^h  )  (2) 


f-iv. 


rhe   s  if::nif  icance 


eqiia  ti  on 


t  i 


O  "; 


is   the   sjistence 


polarisation  independent  ^ulls.  Tne  '^utput  voltif^e  eqijils 
zero  whenever  the  azinuth  an^'le  equals  ?  or  ISC  der^rees  .  The 
incidence  an^le,  the  ■relative  ph^se  and  the  relative 
amplitudes  of  the  electric  fi°ld  conporenrs  do  not  affect 
these  nulls.  It  is  due  to  the  interf erc^^eter  structure  tr.at 
these  nulls  exist;  they  are  therefore  called  interf  ercTie  ter 
nulls  to  distinguish  t  nen  fron  the  si^^ple  loop  nulls.  Ei---;-'-.  re 
3.  tak:en  frorr  reference  11,  ^"aphi  call;.'-  cisplays  tne 
pola  T'i  za  ti  f^n  independent  nulls  f^'"  different  conditions  cf 
i  n  c  i  d  e  n  c  e  a  n  d  p  0 1  a  r  i  z  a  t  i  0  n  . 

Equation  (2^  is  an  expression  for  a  fixed   orientation   of 
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tne  spaced  loop  array  in  the  coordinate  system.  To  rnake  the 
orientation  artitrary  the  variable  alpha  is  introiuced  into 
(2). 

\aO<.  =    -2'^  SIN   SI:\'2(<$)  -  OC:^  +  ^h  3Ii\26  3IN  [  ^  -  (^  (  3  > 

( ^  -  ^"i  :    relative  azinuth  an-.^le 

The  alpha  variable  permits  the  expr-^ssion  of  the  outpi-t 
volta^^e  of  a  spaced  loop  array  oriented  alpha  de^-^rees  frcm 
the  x-ajis  to  be  written  as  shovn  in  \Z)  .  This  v/ill  later 
allow  equation  (3)  to  exp'^ess  the  output  vol  t-3r--;es  of  more 
than  c^.e  pair  ^f  spaced  loops  set  at  different  anrjles  in  the 
coordinate  systen. 
B;'''  definin.^: 


C   = 
A0    ■■ 
A2    = 

T5P       : 


:v   SIN  6 

Eh   SI\(2  6  'exp'  :  ^  h  '  V2 
-C    SI\'   2<^     -   A'^    CCS    2^ 
-A:'    SIN    2<i>     -    C    CCS    2^ 


Equation    (3'    can    be    vritten    as: 


Vaoc     =    A0    +   A2    COS    20<.    +    B2    SIw    2<^ 


(4) 


This  lorr  "oernits  a   Eourier   series   interpretation   cf   tne 
spacea  loop  output  voltd.?e: 

AZ  =  dc  term  of  the  output  volta.^e 
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A 2  and  32  are  coefficients  of  tne  secord  ha^nonic 


The  significance  of  (4'*  is  tnat  fi^r  a  lixed  value  cf 
tarf:et  azir^ith  and  elevation,  the  spaced  loop  voltage  as  a 
function  of  relative  azimuth  is  lirited  to  a  second  narnonic 
of  the  relative  aziruth.  The  application  of  the  Nyquist 
sampling  criterion  ^e veals  that  the  voltage  patte^^n  can  te 
duplicated  "bj  four  sanple  values.  '^herefo^e.  a  spinnir..-; 
spaced  loop  can  "be  synthesized  "by  a  ninimum  of  four  samples 
taKen  equally  spaced  through  35?  de^^rees  of  azirutn. 

The  solution  for  the  cearin.?  (the  azimuth  an.-^le  pnij  is 
derived  from  equation  (4)  and  the  definitions  ^iven  above  for 
C,  A'Z' .  A2  and  '"^Z .  P/  al-=!ecraic  rr.ani  pula  t  i --^n  it  is  letLTr.ined 
that , 


C  =  +/-  (A2  ^  E2  -  A 


(5^ 


^.Z,  A2  and  32  vill  "be  shewn  to  "be  measureacle  c/iantities.  C 
is  determined  frc^  equation  (S^i  above.  Jsirg  tne 
relationships  , 

S  I M  2  4)  =-[!/(  C   +  A ^   )  1  \     (  C  }  ( A2  )  +  ( A^  ;  (  32  )  i     (  o  - 

CCS  2<^  =  -[  l/(  C   +  A.?   )]  f  (A0)(A2;  -  (CM  32  M     '  ?  ' 


one  can  determine  the  azimuth,  phi.  by 


(b  =  (C-.5'  ARC^AM  (  SIM  2  (J)  /  COS  2  (J)   ;   ^  n  •lc2 


n=J,l   V 


Inherent  in  equation  v4:)  and  made  obvious  in  equation   ;3i 
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are  four  null  an^tirriit  ies .  There  are  tvr  nulls  IQ'Z  decrees 
apart  that  can  he  attributed  to  the  si.^-ple  loops  and  two 
nulls  18*3  degrees  apart  that  are  the  interf e'^onieter  rulls. 
The  SV.rl  analysis  shows  that  hy  aid  in:;  tne  sinple  I'^'^p 
phasors  into  the  analysis.  tne  simple  loop  nulls  can  be 
determined  and  then  discarded.  By  conparing  the  spaced  loop 
output  to  the  sinple  loop  cutifut.  the  c^r^ect  interferometer 
null  which  represents  the  desire"*,  bearin^^  can  be  identified. 

The  en.-^ineers  at  S'''^I  used  the  ideas  they  developed  above 
to  design  a  fixed  spaced  loop  array.  Tne  undesirable 
rechanical  feature  of  the  rotatin>=;  spaced  loop  was  elirinateci 
by  usins  four  spaced  loops  fixed  in  an  array  to  synthesize 
rotation  as  shown  below. 


CK 


-t"' 


190 
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Figure  "^ 
Spaced  loop  array  gecmetry 
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The  !\ynijist  critericn  requires  a  rininiun  of  four  sanples  to 
syntnesize  equation  (4).  This  could  reliatly  be  aocomplished 
oy  three  nairs  of  spaced  loops,  tut  to  provide  for  additional 
reliatility  in  the  presence  of  noise,  a  four  pai:-  spaced  loop 
array  was  constructed. 

Assur^i^g  the  orientation  ^iven  in  the  dia^ran  aVove,  one 
can  deter-^ine  A0,  A2  and  "32  in  ter-ns  of  tne  iniiviiudl  spaced 
locus  .  Solvin:^ 

Vao^  =  A'?  +  A2  CCS  2  oc  +  32  SIN  2  (X 

in  terns  of  al-oha  yields, 


(?(s      = 

2 

o(.  = 

45 

oc  = 

92 

0^  = 

13 

Vai)  =  A^  +  A2 
Va45  =  X2    +32 
^1^92    =  A^  -  A2 
Val35  =  A0  -  32 


where  Vav" .  Va45 ,  Va90  and  Val35  are  the  Phasors  of  tne  spaced 
loons  in  the  array  shown  above. 

Tnis  provides  four  equations  to  solve  for  three  unknowns.  One 
solves  for  A0.  A2  and  32  oy  the  following  equations, 

A0  =  (2.25)     (Vae  +  Va45  ^    \la9d    ^    7a  135' 
A2  =  ^^.5)  (Va?  -  Va90) 
32  =  (^.5)  (Va45  -  7al35) 


After  A?,  A2  a^-d  32  a'^e  deterrnined  fror^  the  phasor   equations 
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atove,  they  are  substituted  into  equations  (5'  Ar.±  (7;  vni  oh 
in  turn  are  used  to  solve  equation  (8)  for  the  i our  possible 
bearings . 

?urtner  al.^ebraic  and  tri-^oncnet  ri  c  analysis  detailed  in 
Reference  [9]  shows  tnat  the  sinple  loop  nulls  can  te 
de terriined  by , 

CX  -  AP.CTAN  (  VL?/-VL92  ^  =  cj)  -  A^.CTAN  (  C/.i^  ) 

vnere  VLZ  and  VL9'2*  are  the  Dhascrs  cf  the  two  simple  loop 
Dairs  . 

Once  the  sinple  loo"o  nulls  are  .-en own  .  tne  si^n  of  C  in 
equation  (5''  can  be  dete"^nined.  This  in  turn  leads  to  ihe 
unanbi^'uous  selection  of  the  nron>er  interf erof^ieter  null. 


4)  =  APCTAM  (  VI /Y2 


wnere. 


VI  =  [j/'f-A:^  -CM  f(C)(VLP^)  -  (A5)iVLy? 


V2  =  f.1/(-A0^-  C*^)]  \     (C)(VL3)  ^  (A0)(VL9<3) 


It  was  noted  above  that  four  equations  are  available  to 
solve  fo"  three  un'-mov.'n  coefficients.  Tne  additional 
information  pe^r^its  two  separate  solutions  fc  the  A'<?  tern. 


fA^l  =  (7.5'     (  Va?  ^  '/age  '• 

fA'Z]  '  =  f  7.5  ■  (  Va4:5  -^  "7al35  ) 

The  '^.i  f  f  e""ence  tetween  these  two  ^.^  ter'^s  should  ideallv  be 
zero.  If  the  liffe^*=nc:e  is  rot  zero  the'^e  is  an  inconsistency 
within  the  syster".  This  difference  is  calle'i  tne  A4:  tern 
"because  it  corresponds  exactly  with  the  coefficient  of  the 
fourth  har'-onic  of  a  Fourier  series  expansion  of  tne  spaced 
loop  pattern  in  azir^uth.  The  A4:  ter-n  is  therefore  en 
important  para^ieter  in  determining  oearin.^  viuality. 

C.  SYSTEM  DiSl^N  AMD  IMS'"F.UMSNTATIOi\ 

The   aT'ray   of   spaced  loops  and  simple  loops  suitaole  for 
past  mounting  shown  i^  fi5:ure  8  was  built  cy  S.vRI.  There   are 


four   pai^s 


spaced   loops   in   the   lover  bay.  ?^ach  pair 


consists  of  4'^  inch  high  tv  2?  inch  wide  simple  lo^ps 
separated  by  60  inches.  The  output  of  tnese  diamet  T-ical  ly 
opposite  loops  are  connected  in  paT'allel  opposition.  The 
simple  loops  in  the  upper  bay  are  used  to  '^esolve  ambiguities 
in  the  bearin.^  alf^orithm.  These  diametrically  opposite  simple 
loops  are  connected  ii  parallel  assistance.  The  reference 
antenna  is  synthesized  by  quadrature  addition  of  tne  simple 
loops . 

A   block  diagram  of  the  equipment  suite  is  drawn  in  fi^^ure 
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Figure  8 
Coaxial  Spaced  Loop  Mast-Mountable  Array 
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Figure  9 
Block  Diagram  of  SWRI  Spaced  Loop  Antenna 
System  Instrumentation 
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9.  The  RJ  seqvercer  is  a  conpiiter  controlled  switch  that  is 
necessary  to  provide  oiigh  speed  sv/itching  of  the  different 
elenents  of  tie  array.  A  dral  channel  receiver  is  used  to 
r)rovide  3  receiver  channel  for  the  reference  si.^nai  and  a 
channel  for  tne  loop  volta^^es.  Tne  predetected  output  of  botn 
channels  of  the  receiver  is  monitored  by  a  mnase'neter  ihat 
provides  a  di2:it3l  ^easurenent  of  tne  n-ase  of  the  antenna 
elements  with  respect  to  the  ^efe^ence.  The  detectors  are  a 
pair  of  p-^ecision  pea'^  detectc^s.  The  detector  output  is 
sappled  and  di.-?itiz°d  hv  the  analo.-^  to  di^jital  con  verier.  The 
digital  data  is  ^oute^^  to  a  niniconputer  for  process  int<: .  iacn 
data  frar^e  is  apprcxinatelv  2'?  "ns  in  duration.  Tne  data  frane 
consists  of  tie  six  co ropier  numbers  representin-^  tne  six 
voltage  phaso^s  'Va-^.  Va4:5,  Va90.  yal35,  VLZ  and  VL9^  '  .  Net 
all  of  tne  data  frares  a^e  acceptable.  Tne  vcltage  values 
nust  be  vithin  the  linea"  "^an^'^e  of  the  detector  and  receiver 
circuits.  The  acceptable  <i&t^  franes  s'^e  tne  inp^:t  to  tne 
algorithms  that  detp^nine  the  beaT'inrs  an'"  resolve  tne 
ambiguities  . 
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IV.   DATA 


A.  DATA  FILES 

The  SW^I  equipment  suite  is  aTran^ed  so  that  data  measured 
from  the  spaced  loop  antenna  instrumentation  is  stored  on 
magnetic  disk.  This  permits  the  DF  operator  to  postprocess 
the  data  usin^  statistical  tecnniques  to  derive  a  more 
accurate  tearing.  Also  availalle  is  the  capatility  of  mass 
storage  on  magnetic  tape.  It  vas  on  magnetic  tape  tnat  SV.RI 
provided  the  Naval  Postgraduate  Soncol  -.vith  nine  files  of 
data  in  1979  and  four  files  in  I960.  The  1979  data  consists 
of  the  fcllcving  files.  The  test  source  vas  a  transmitter 
placed  close  to  the  array  to  provide  a  ground  wave  in 
approximately  the  same  direction  as  'V.vV. 


File    Num 

he  r 

Free ue 
('MHz  ) 

TiCJ 

Source 

Time 
(CST) 

Bate 

1 

10 

l^WV 

28:30 

2/13/79 

2 

20 

VWV 

12:15 

2/13/79 

3 

15 

>/'aV 

I2:5d 

2/13/79 

4 

5 

wwv 

?7:?0 

2/14/79 

5 

15 

wwv 

09:o0 

2/l±/  79 

e 

5.(' 

'1 

Test 

7 

10. 

01 

Test 

e 

15. 

01 

Test 

9 

2d. 

01 

Test 
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The  1983  files  are  as  follows: 


File  \'urber  Frequency 
(MHz  ) 


S  0  u  r :;  e 


•T-i 


T  rni 


Date 


(  CST  ) 


1? 
11 
12 
^7^ 


5 

1^5 

15 

t  .666 


u  H  Y 
WIVV 

'/;wv 

KLC 


22r40 

07:4:0 
09:20 
0£  :15 


2/6/5^ 
2/5/30 
2/5/c'? 
2/7/80 


It  is  known  that  file  1:^  is  a  low  3NR  data  set.  Files  11   and 
12  a  "e  data  sets  with  S\'H's  in  excess  of  2?,   dS. 

B.  'a  Pi    AND  KLC 

'aV'V  is  an  ideal  target  tecause  it  is  an  anplituie 
nodulated  sip^nal  with  no  carrier  suppression.  For  the 
majority  of  the  hourly  information  duty  cycle,  tne 
information  modulated  on  the  carrier  is  simple  4 4: J,  ^ii  and 
600  H?  tones;  ticks  I  and  occassional  voice  announcements.  The 
WWV  sif3;nal  is  stable  to  +/-  two  T)arts  in  \  1<  ;  .  and  it  is 
availacle  on  5  ,  10,  15  and  20  MHz  for  24  dours  a  day.  Tne  v/',V'/ 
signal  is  transmitted  fron,  Eoulder,  ColoT-ado,  which  is 
geographically  fiied  at  43.9  M  and  105.1  a.  The  true  Gearing 
of  the  £-reat  circle  arc  passing  througn  San  Antcnio.  Texas, 
and  Boulder  is  336.7  degrees,  and  the  length  of  tie  arc  is 
1367  km. 

KLC  is  a  manual  morse  ship-t  c-snrre  station  t  ransmi  t  ti '^ -^^ 
f'-om  a  platform  in  ^he  Grilf  of  Mexico.  This  signal  was  cncsen 
as  a  target  "because  of  the  on-off  leying  (00?I)  modulation  and 
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because  its  relatively  short  ii  stance  fror,  San  Antonio 
results  in  a  s-rywave  vith  a  high  angle  of  arrival  at  tne  S'a'PI 
antenna  array.  The  COK  nodMla ti on  is  important  cecause  it  is 
a  favorite  node  of  tactical  cormuni  cat  i  ons  »  it  is  "rrief  and 
reliable.  Tie  true  bearing  f^orr  San  Antonio  to  r'LC  is  '?S9  9 
degrees  ard  the  distance  is  372  '<.r . 

C.  DATA  RECORDS 

Each  file  consists  of  13,<3^0  records,  fiacn  record  consists 
of  the  A?,  phase.  A4:  and  tearing  terrns  calculated  from  frarnes 
of  data  (six  voltage  phasors)  vnich  were  generated  every  i-j 
ms.  The  A'5  in  the  data  is  a  nornalized  version  of  tne  Ad 
explained  in  the  previous  section.  In  that  section, 

A'J  =  (2.5'  Zh(SIN2  9  )exp(  j  <^h) 

The  data  on  the  tape  is  A?  nor-^ali^ed  by  the  factor  ZvSI.^  . 
which  yie Ids : 

A0n  =  [  Sh/E^,^  1  CCSO  exp(j<i>h^ 

This  is  a  complex  nurnter  of  whicn  only  the  Pagniturie  is  use"". 
All  further  references  to  A?  in  the  data  and  analysis  section 
vill  be  the  r^a^nitude  of  the  nornalized  value: 

A0  =  r  Eh/Ev  1  CCS  6 
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It  should  t.p  n^ted  he^^e  tnat  the  A'3  terr  is  a  neasijre  of  t.ie 
arrount  cf  horizontal  p-^l  a^^i  za  ti  on  present.  If  E.-i  is  greater 
than  Sv,  the  "atin  ^h/?v  will  tend  to  ma':ce  Xi  a  auncer 
greater  than  one.  If  the  anpl'=  f^on  the  perpen-l  icul  ar  ,  tneta. 
is  large,  there  is  a  greater  effective  a^-r-ay  aperture  in  tne 
plane  of  Iv ,  and  A^  is  srall^^.  If  Kd  is  large.  tne 
horizontal  co:nponent  is  dominant.  If  A<?  is  small,  tne 
vertical  component  dominates. 

The  second  term  of  tne  record  is  the  phase.  It  is  tne 
phase  of  the  horizontal  electric  field  relative  to  the 
vertical  component  of  tne  field.  It  is  a  calculated  value 
between  -1-?  and  +1SZ  degrees.  If  this  :)hase  angle  is  a 
constant  ^erc.  the  polarization  is  linear.  If  it  is  a  nonzero 
constant,  the  TD'-la^i  zat  ion  is  elliptical.  ?he  phase  value  is 
typically  noted  to  vary  randomly  within  a  limited  range  OYP.r 
short  tine  durations.  Cver  durations  of  several  minutes,  ii 
v;ill  varv  over  the  entire  -180  to  -^18/"?  degree  range  due  to 
changes  in  Ih   and  Ev  nath  lengths  and  mr^ltimode  interference. 

The  third  term  in  a  record  is  tne  amplitude  of  the  M  term 
which  v/as  discussed  in  the  previous  section.  It  is  a  measure 
of  the  inconsistency  within  the  spaced  loop  H?j?  system.  It 
is  in  large  part  due  to  noise,  but  it  can  also  be  tc  a 
limited  de?ree  a  measure  of  circuit  imbalance,  measureme^'t 
e^ro^.  component  failure,  software  failure  and  site  error. 
Its  value  is  that  i^  is  a  measure  ^f  "oerf 0"mance ;  nowever.  it 
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is  not  a  syster^  diarnostic  tool. 

The  fourth  nenber  cf  the  record  is  the  calculated  bearing. 
It  is  an  integer  valre  fT'cn  ?  to  359  desjrees  .  ?or  the  .v.vV 
signals  this  bea^irs:  is  the  syster;;'s  estinate  cf  the  an.^le  cf 
arrival  of  the  signal  wavefront  which  should  not  vary  f^-^ 
fror  the  value  336.7  degT-ees.  ^or  5CLC  the  true  bearinr  is 
7:89  .9  degT-ees  . 

D.  ICNO SPHERIC  DATA 

N'o  ioncsTJheric  sounding  infor~iation  v/as  available  for  tne 
tine  periods  during  which  tne  data  was  recorded,  iowever, 
pronagaticn  information  was  nrovided  by  tae  Naval  Gcean 
System  Center  (NOSC''.  Using  tne  Known  sun  spot  nunber  for  tn.e 
data  recording  dates.  they  employed  a  computer  propagation 
prediction  pr'^gT-am  known  as  PROP'-I'^T  to  Droviae  ray  trace 
diagrams  J  MU7 .  LU?  and  powe"  predictions!  and  24  hour  line  cf 
bearing  variance  cu^vps  .  This  data  represents  a  good  estimate 
of  propagation  conditons  between  Boulder  and  San  Antonio  f;r 
the  times  of  interest.  Examples  of  the  program  output 
nrcaded  bv  N'OSC  are  reviewed  in  Appendix  C.  L'sing  tne  NCSC 
data,  propagation  information  for  the  \ 'iV  files  is  taculated 
in  TABLE  I . 
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TABLE    r 

3iin?marv   of    I  onos^'ieric  Dai  a 

File  ?req   Ti^e    Date    Variance      MU7      LU?  Ionospheric 

if  MHz      GMT  d.e,^rees        MHz      MHz  Mode 

5  a  u  a  r  e  d  ?  r  e  a  i  c  t  i  o  n 


1  1?       14  •3'^    2/13/79      1  17 


2  20       ie:15    2/13/79      1  22        5 


3     15   1£:3'J  2/13/79   1      Z2 


13-00  2/14:  ^79   3      12    2 


15   17:3?  2/14/79   1 


I'd        5    2:1(7   P/*-*/?'. 


9P     c^ 


12 


11    ID   13:1:?  2/5 /£^    3      1?    2 


12    15   17:0^  2/5/92    1      23    5 


2.5  Drctable  multimode 
1,  2,  ^3  aops 
possible  multinode 
1  hc^ 

hif?il7  pro  cable 
nultimode  1£.2  nops 
possible  ^ulti'ncde 
ter'~'inatcr  1.2.3  hops 
probable  single 
r.  0  d  e 

possible  rultinoda 
1  .26^3  hops 
P'-n table  mill  t  imo  ie 
terr^inato^  1,2  hops 
sir.csle  noie 
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V.  A^'ALT5IS  OF  SWPI  DATA 


A.  INTRODUCTION 


The  purpose  of  this  analysis  is  to  erar  tne  capafcilitv  of 
the  S'"^?I  spacel  Iood  h;?j~  antenna  s'^sten  to  ieternine  tne 
angle  of  arrival  of  a  short  dvration  sif^-nal.  A  short  duration 
signal  is  consiierei.  to  be  fro-n  r?<J  to  l^dd  "'s  in  dn-'dtion. 
It  is  imp'^rtant  to  note  that  this  is  not  a  f^eneral  analysis 
of  the  perforriance  of  the  antenna  system,  tne  bearinf^  a.id 
sense  algorithrs  or  the  post-processing  al.gor i tn.^s  developed 
by  SV/FI.  It  is  also  inn)ortant  to  recognize  tnat  s/sterr 
develcDrent  is  not  conplete.  but  is  tne  subject  of  on-goin>j 
research.  Tne  data  provided  to  this  investigator  'Aas  provided 
from  a  system  configuration,  not  optimized  for  short  durati:tn 
signals  or  for  sore  of  the  tar_?et  frequencies  reccdec^  . 
Because  the  subject  of  this  analysis  is  short  d'lration 
Si<:yv£ve  signals,  the  following  analysis  takes  into 
consideration  the  need  to  naive  raxinur:  use  of  tne  available 
riata.  V.' he  re  as  3V;P.I  algorithms  stress  a  bearing  selection 
process  that  eliminates  a  large  percentage  of  trie  data 
records  tc  enhance  the  reliability  of  the  esti-atec  bearing, 
this  analysis  reccrnizes  that  a  si^:!nal  of  £00  ms  duration  is 
represented  bv  only  1'?  data  frames  and  that  some  compromise 
to   reliability  must   be   made.  Tne  term  reliability  in  tnis 
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report  is  used  as  a  neasn^e  of  confiience  in  the  validity  '^f 
the  data.  If  one  associates  a  standa-^d  deviation  of  53 
de-^rees  with  a  data  record  and  2?  de^^rees  vita  a  second  data 
reco^'d,  the  latter  would  be  considered  rc^e  reliable. 

The  analysis  fi^st  concentrated  on  examinin:^  the  data 
pnirna^ily  h''  filte^ins*  on  the  A 4  te^n.  the  indicator  r? 
syster  inconsis  te^^. c''.  If  the  Ai  terp  is  s'^all,  the  bearing  in 
that  record  should  "be  considered  more  reliable  tnan  a  oearinv^ 
associated  with  a  high  A4  value.  Using  tnis  approach  a 
fCRTPAN  pro>=!rani  called  D??PP  ( D?  E^?or'  was  developed  to 
examine  each  file  and  report  the  average  hearing  error, 
standard  deviation  of  bearing  error  and  two  otner  statistics 
of  short  duration  sifrnals  with  respect  to  tne  Ai  term.  A  full 
discussion  of  the  DT'ogram  is  detailed  in  the  next  section.  It 
was  discovered  that  tne  A-l  te  r-n  is  a  useful  parameter  for 
determining  hearing  reliability  in  tne  majority  of  cases. 
However,  when  the  A^  th'^eshold  is  set  to  only  allow  the  data 
record  with  A4:  approaching  close  to  zero  (the  tneoretical 
ideal).  there  is  not  a  large  probability  -^f  determining  a 
bearing  on  a  short  duration  signal.  In  an  attempt  to  improve 
on  the  sole  use  of  the  A4  term  as  a  reliabilty  inaicator,  a 
probabilistic  likelihood  ratio  matrix  based  on  all  of  tne 
available  signal  parameters  was  employed  in  the  analysis.  Ihe 
details  of  this  approach  are  in  the  section  titled  LriAT 
(Likelinood  'iATrix). 
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A  closer  exan^i  nei  na  ti  on  --/as  '^aae  of  the  Drctle-n  of  bearir^^ 
ambisiuitv.  The  tech'^. ique  employed  in  this  po"tion  cf  t.ie 
analysis  is  ^^iven  in  the  section  titled  A^-liilGUIT'f  ?. l^iSOLUTICN . 
This  analysis  r;i"'^es  some  useful  insi,^ht  into  possible 
difficulties  within  the  =ntenna  systen  that  may  prove  to  te 
the  most  tracts  Die. 

B  .  DFERR 

The  purpose  cf  this  analysis  vas  to  determine  nrv 
accurately  a  DF  tearing  could  be  calculated  f'^or  the  :^iven 
data.  The  data  consists  of  230  seconds  of  WaV  per  file  (one 
file  of  KLC).  To  study  short  duration  signals  it  is  only 
necessary  to  consider  the  2d'Z  seconds  of  data  to  be  a 
continuous  concatenation  of  short  duration  si.-^nals.  Tne 
l'd,2Z'^  records  in  each  file  contain  tne  data  for  t.ie  2'^  ms 
sa^iplin^  T)eriods»  therefore,  integer  multiples  of  records 
correspond  to  different  si-::nal  duratio^^s.  To  eramine  systen 
performance  against  a  2<.V  ms  signal,  one  need  only  examine  a 
file  1?  records  at  a  time.  A  20c^  second  file  may  be  thcu.g'.it 
of  as  containins^  1?33  si.^nals  of  23^  ms  dur^ation.  Similarly, 
for  a  signal  duration  of  1  second.  53  records  rr:ay  be  used  to 
synthesize  the  signal,  and  the  file  is  made  up  of  2c'^ 
s  ignals . 

A   FORTRAN   tro^^a^  named  DFFRR  v/as  v/ritten  to  examine  tne 
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data  files  based  on  the  above  conceT>t  .  The  prO;?ran  vas 
desii^Tned  to  examine  sip-nal  durations  frorr  20  ins  to  200 
seconds;  however,  it  was  used  for  this  analysis  in  two 
ran^^es.  l^'d  ns  to  1'?.'0  "^s  in  increrents  of  12"?  rs  and  J. 
second  to  1?  si^conds  in  increments  of  1  second. 

The  fireneral  purpose  of  DFTFF  is  to  examine  system 
performance  as  two  parameters  are  varied.  The  first  parameter 
is  signal  duration?  the  second  is  the  A't  term.  Tne  Ai  term, 
explained  in  section  III,  is  the  measure  of  inconsistency 
within  the  D?  svst»m.  If  ^he  Ai  term  is  lar^e.  tne  bearin,^ 
value  in  a  "ecord  is  not  considered  reliable.  If  aI  is 
small. more  confidence  is  "olaced  in  the  bearing.  The  relative 
descriptors  larsre  ^v.''.  small  have  vet  to  be  evaluated.  In 
order  to  evaluate  the  pertinent  "an>^e  of  k^  values,  the  A4: 
threshold  'A4:>'AX)  is  varied  between  a  small  value,  ?• .  1  ,  and  a 
larse  value,  1.2,  in  increments  of  .! .  1 . 

If  the  value  of  A4:  in  a  record  is  ecuel  to  or  less  tnan 
the  value  of  A'4:MAX  set  in  the  pro^;ram,  tne  bearinr;  is 
consi'^.ered  acceptable  and  used  in  further  statistical 
iDrcces  5in^^.  If  tne  value  of  Pi^  is  above  the  limit,  the 
bearing  of  that  record  is  discarded. 

An  explanation  of  fu^tner  DFERR  processing  is  best 
presented  using  an  example.  Suopose  that  the  following 
records  are  being  urocessed. 
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hecord  ft 

A '3 

Phase 
(de-) 

A  4: 

Bearing 
(dec^ 

•  •  • 

•  •  •  • 

•  •  • 

350 

^.821 

-47 

J.  132 

333 

351 

0.611 

-«0 

3  .315 

540 

352 

0.432 

-55 

3.115 

^<lo 

353 

e.512 

120 

3.413 

253 

354 

0.315 

-80 

3.176 

233 

355 

1.311 

-23 

3.215 

343 

•  -  • 

.  •  •  • 

.  .  .  •  . 

•  *  • 

The  value  foT-  the  signal  d^^atio^  is  133  rs  arc  the  A4:'^'AX 
value  is  3.21  the^efo^e,  re::o"ds  353  th^\i  354  are  examined  as 
representing:  a  signal  of  133  ms  duT'ation.  F.ecoT^d  355  is 
irr,media tel.7  rejected  lecause  the  value  of  A4  is  r:reater  tnan 
A4MAX.  The  remainins:  fou^  ''eco^ds  a^e  oalled  "Ai  adnissitle" 
and  are  used  to  determine  a  tearing  nean  and  standaT-d 
deviation:  (The  rrean  and  standard  deviation  formulas  used  in 
DFERR  are  derived  in  A^Deniiy  A.) 


MTAN  =  3  3?  dec 


3TD  =  44  de;^ 


These  statistics  ar°  used  to  foT-m  a  v/indow  centered  at  339 
degrees  extending:  44  decrees  on  eithe^'  side  of  the  mean.  The 
A4  admissiole  hearings  must  pass  through  this  window  lo.r' 
further  consideration.  Pecord  numher  351  is  not  '  windcv; 
admissitle"  an'"*,  is  discarded.  Th^'  remainin^T  tr.ree  records, 
"bein^  hcth  A 4  and  windov/  admissible,  a  ""^^  used  to  compute  a 
second  mean  ani  standard  deviation: 


MFAM  =  335  '^eg 


3TD  -  4  leg 
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This  rean  is  DFilRH's  test  estinate  cf  tiie  bearing?  for  this 
one  1"?^  PS  si-?nal  f^-ecords  ZbZ    throi^^^'f^  354  •  . 

This  reported  bearing  is  con-pai^el  with  the  tri:e  Dearin^;. 
3.37  'ie^,  an'.''  the  hearing  error  is  ^onpi^ted  as  335-337=-2  deg. 
Additionally,  the  valid  sir^al  counter  is  augmented  by  one. 
The  nurber  of  valid  signals  will  be  used  later  to  deiernine 
the  probability  of  obtaining  a  tearin<?:  (?0:). 

If  there  are  less  than  three  A4:  adnissible  records  or  less 
than  two  window  adnissible  records  in  a  given  signal 
duration,  the  signal  is  considered  invalid  di.e  to 
insufficient  '^ata  and  is  counted  in  an  invalid  signal 
counter.  To  avoid  the  loss  of  reliable  data  in  the  case  of  a 
snail  standard  deviation  of  the  A4:  adnissible  beamings,  the 
screening  window  is  not  permitted  to  be  narrower  than  Ic 
degrees . 

The  D?5;?.R  Drogran  processes  ten  separate  signal  durations 
at  ten  aifferent  A4MAX  settings.  Pro-ran  output  consists  cf 
four  tables  on  separate  "oa^-es  .  tach  table  is  a  ten  by  ten 
matrijJ  the  rows  correspond  to  tne  A4MAX  vlaues  and  tae 
columns  correspond  to  the  signal  duration,  see  tables  II 
through  7.  The  first  table  (table  II)  is  the  average  bearing 
error.  This  is  the  average  of  all  the  separate  means 
reported.  In  the  example  above  for  if^o  ms  ,  the  mean  335  would 
be  one  of  a  maximum  possible  2/;2!?^  values  that  would  be 
averaged  and  then  re"Dorted   in   tae   row  rA"Aky.-c.l      and   the 
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colnnn  sigral  du^'ation  equals  IP'^  ns  . 

The  second  tatle  (table  III)  is  the  standard  deviation  of 
the  bearing  e^^ors  reported  in  the  first  tatle.  ^or  tne  l'<^'i 
ms  signal  duration  category  the  standard  deviation  vould  re 
of  a  maxirum  possible  23^0  mean  values  (the  actual  sample 
size  is  ecual  to  the  number  of  valid  sirnals).  As  tne  si;::ndl 
duration  increases,  the  sample  size  decreases.  For  tne  lJc:0 
ms  column,  the  maximum  possible  number  of  si.-;nals  is  22^.  If 
medium  duration  signals  are  examined  with  D^iRR,  one  must  be 
attentive  to  the  sample  size.  For  a  Id  second  duration  si.^nal 
cate^rorv,  the  samule  size  has  diminished  to  twenty .  For  l'/-^ 
second  signal  durations,  there  are  only  two  samples  and  tne 
validity  of  a  standard  deviation  is  highly  questionable. 

The  third  table  of  output  (table  17)  is  tne  avera,f.;e 
intra-signal  standard  deviation.  This  is  tne  average  of  tne 
standard  deviations  reported  for  individual  signals.  In  the 
previous  example,  the  5TD=4:  vould  be  one  of  a  raximum 
possible  260?  standard  deviations  to  be  averaged.  If  tne 
average  intra-signal  standard  deviation  were  equal  to  four, 
the  interpretation  would  be  that  for  all  signals  of  I00  ms 
duration,  after  the  unreliable  bearings  are  discarded.  the 
expected  standard  deviation  of  the  remaining  cluster  cf 
bearings  is  four  degrees.  The  term  intra-signal  is  used  to 
distinguish  it  from  the  standard  deviation  of  bearing  errcs. 

The   fourth   table   (table   V)  is  the  compilation  of  valid 
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signals  for  each  of  the  A4:MAX  and  sir^nal  duration  categories. 
This  is  an  important  statistic  needed  to  compute  the 
protatility  of  obtaining  a  bearing.  In  the  case  of  tne  I'^d  rs 
signal  duration,  2000  are  possible.  If  A4MAX=t:/.2  and  12e'0 
signals  are  valid,  tne  FOB  for  l<d^  ms  is  equal  to  the  ratio 
of  valid  signals  to  possible  signals.  In  this  case.  F0B=3.6. 
In  general,  as  either  A-i^IAX  or  tae  signal  duration  increases, 
the  number  of  valid  signals  approaches  the  number  of  possible 
signals.  For  A-^MAX  above  3.4-  cr  the  si.^nal  duration  above  one 
second,  the  POE  is  approximately  one. 

DFEFP  was  used  to  process  all  of  the  files,  including  tne 
OCK  modulated  KLG  signal.  The  rLC  file  required  a  slightlv 
modified  version  of  DF^RH  because  each  record  of  data  was 
searched  fo^*  a  flag  that  indicates  that  the  signal  is  indeed 
present.  SeaT'ch  was  also  inaie  fo^  flags  that  indicate  that 
tne  system  is  saturated.  KLC  is  considered  a  file  of  signals 
separated  bv  noise.  It  is  not  a  continuous  signal  lii^e  the 
WWY  signal;  in  fact,  the  duty  cycle  is  less  than  twenty-five 
percent . 

The  tabular  output  produced  by  DF5RR  can  be  considered  a 
second  level  data  base.  To  simDlify  follow-on  discussions, 
the  statistics  in  the  DF2RR  output  will  be  referred  to  as  tne 
"four"  statistics.  5:xamination  of  the  DFi:RR  data  base 
revealed  that  the  FOU^  statistics  are  a  strong  function  of 
the  A4  term.  As  the  A4  term  increases,  corresponding   to   tne 
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acceptance  of  more  inconsistent  data,  the  standard  deviation. 
intra-si,^nal  standard  deviation  and  nunter  of  valid  signals 
increases.  Based  on  the  lar^e  standard  deviations  recorded 
above  A4=C4:  and  the  insignificant  increase  in  PC£  atove 
A4=«2.4,  the  range  of  interest  was  restricted  to  A4  values  in 
the  0.1  to  ^^.4  region.  Even  witn  this  restriction,  tne  arrount 
of  data  is  tco  lar<?e  to  present  in  this  repcrt;  however,  it 
should  he  noted  that  the  changes  in  the  FOUR  statistics  ere 
typically  monotonic  as  A4  varies.  The  data  "oresented  herein 
is  for  A4MAX=0.2  and  A4MAX=e.4.  The  best  Derfcrnance 
category,  A4MAX=0.1.  is  not  Dresented  because  of  its  low  PCS 
and  because  it  is  not  significantly  different  from  A4i'lAX  =  'i? .  2 

The  tabular  data  of  DFSRR  does  not  permit  easy  visual 
perception  of  the  characteristics  of  the  FOUR  statistics.  A 
plot  prosran  was  written  to  display  the  data.  Tne  figures  at 
the  end  of  this  section  are  of  the  wrvV  5,1?  and  2v?  MHz  files 
of  19'^9  (Fig.  l'?-24)  and  198v^  (Fig.  25-3'i).  For  eaca 
frequency  the^e  ave  five  graphs.  Ir.  eacn  set  of  five,  tne 
first  two  and  the  last  are  the  most  important  as  tney  are 
concerned  with  the  short  duration  sirnals.  The  second  two 
graphs  are  fo^^  signal  duratiors  of  one  to  ten  seconds.  They 
are  included  to  illustrate  how  the  FOUR  statistics  tend  to 
reach  steady  state  values.  The  fifth  graph  is  the  histogram 
of  the  bearing  errors  for  J\4MAX=0.2  and  signal  duration  equal 
to  22d     ms .   Annotated  in  the  upper  right  hand  corner  of  all 
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the  corrpute^  drawn  ^raDhs  are  the  average  bearing  (3  . 
standard  deviation  (-')  and  avera.^e  bearin^^  error  (  ~J  for  the 
entire  file  portrayed  b?/-  the  eraph.  These  are  tne  statistics 
foT  a  si£!nal  duration  of  200  seconds  subject  to  the  A4>1AX 
constraint  labelled  on  the  e;raph.  If  one  assunes  that  tne 
true  angle  of  arrival  was  fixed,  that  no  multipatn  existed 
and  that  there  was  no  slow  tern  variance  due  lo  tne 
ionosphere.  the  curves  plotted  should  converge  to  tne  2<''>i 
second  average  bearing  error.  Both  the  standard  deviation  and 
the  average  intra-signal  standard  deviation  should  converge 
to  the  20?  second  standard  deviation.  The  P03  should  converge 
to  one.  However,  the  PROPHET  program  graphs  show  that  in  many 
cases  there  probably  existed  multimode  conditions  (actual 
amount  of  interference  un'tcnown )  and  that.  in  all  cases, 
varying  degrees  of  variance  existed  due  to  polarization 
fading  and  interference.  Despite  this  difficulty.  the  220 
second  statistics  car  be  considered  an  approximate 
convergence  point. 

The  first  curve  (X)  of  each  graph  is  the  average  bearing 
error.  This  is  the  difference  between  the  calculated  bearing 
and  the  true  geographic  bearing  to  the  signal  transmitter 
(337  or  090).  The  most  prominent  feature  of  all  the  average 
bearing  error  curves  is  that  tnev  are  not  a  strong  function 
of  signal  duration.  This  could  aave  been  anticipated 
realizing  tnat  the  average  of  the  majority  of   small   subseis 
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of  a  lar^e  set  will  closely  approximate  the  avera^re  of  tne 
large  set.  In  this  respect  the  avera.^e  tearing  error  as  a 
function  of  signal  duration  is  not  particularly  useful.  A 
much  more  significant  view  of  the  average  bearing  error  is  a 
histogram  of  the  hearing  errors  from  each  valid  signal.  A 
histogram  for  A4MAX=0.2  and  a  signal  duration  equal  to  2Z0  ms 
is  the  fifta  figure  in  each  set  of  five  ^Fig. 
14, 19,24,29,34) .  These  histograms  show  that  the  distritution 
of  the  bearing  error  for  mcst  of  the  files  is  only  rougnly 
approximate  to  a  normal  distribution.  I^ajor  deviations  from 
the  normal  curve  are  the  accumulation  of  bearings  in  tne  leO 
degree  ambiguity  rerAon  and  the  dis  t ''ibuticns  wnere  multimcde 
propagation  was  highly  probable.  The  major  difficulty  with 
the  average  bearing  e^ror  data  is  the  absence  of  calibration 
data.  Tae  amount  of  correctable  bias  error  is  unknown. 

The  standard  deviation  curve  (''M  is  tne  most  important 
informaticn  on  the  graphs.  If  a  bearing  is  to  be  used  with 
other  bearings  to  compute  a  fix,  the  standard  deviation  is 
used  to  compute  the  fix  area  for  a  given  probability  tnat  the 
target  will  be  within  the  fix  area.  The  standard  deviation  is 
also  a  measure  of  confidence  in  a  single  bearing.  If  the 
bearing  distribution  is  normal  and  the  standard  deviation  is 
2z  degrees,  one  can  expect  that  the  bearing  calculated  is 
within  2'<?  degrees  of  the  true  bearing  (assuming  aveT'age 
bearing=true  bearing)  about  67  percent  of  the   time.   Because 
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the  distributions  for  the  data  files  are  only  roughly  normal. 
amplification  of  the  significance  of  the  standard  deviations 
is  pertinent.  The  following  data  are  the  approximate 
percentages  (within  five  percent)  of  bearings  falling  witnin 
the  standard  deviation  of  each  file  for  a  2^'i  ms  signal 
duration  and  AiMAX=0.2. 


File  # 


1  (Fig.  19^ 

2  (Fig.  24) 
3 

4  (Fig.  14) 
5 

10  ^Fig.29) 

11  (Fig.  34' 

12  (Fig.  35^ 
13 


STD  % 

31.1  90 

25.6  35 

31.1  75 

15.5  55 

37.4  90 

59.6  75 

29.5  95 

16.2  39 

19.3  33 


This  data  indicates  that  the  distributions  are  denser  than 
the  normal  and  that  system  performance  is  better  than  one. 
thinking  in  terms  of  a  normal  distribution,  would  believe. 
For  example,  in  file  11  (Fig.  34)  the  standard  deviation  is 
about  30  degrees,  but  95  percent  of  tne  data  are  within  tne 
standard  deviation.  The  expression  of  standard  deviation 
cannot  be  separated  from  its  distribution  and  still  retain 
meaning . 

Observins:  the  curves  [Yie.  10-13,  15-18,  -c0-23,  25-28, 
30-33^  the  reader  will  note  that  the  standard  deviation  is  in 
almost  all  cases  a  monotonically  decreasing  function  of 
signal  duration.  An  interesting  comparison  is  the  standard 
deviation  at   100  ms  (STDl)  and  at  10  seconds  (STD2),  again 
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witn  A4MAX=?.2 


Pile  ^   3TD1   STD2   3TDl/srD2 


l" 

29 

5 

5.3 

2 

25 

11 

2.3 

3 

31 

13 

2.i 

4. 

18 

5 

3.6 

5 

42 

10 

4.2 

13 

57 

16 

3.2 

11 

30 

4 

7.5 

12 

16 

9 

1.8 

13 

17 

6 

2.3 

These  results  are  conparatle  to  those  reported  by  Eiin 
[Ref.3l  and  consistent  in  order  of  magnitude  to  the  tirrie 
required  to  average  the  fluctuating  surfaces  of  constant 
amplitude . 

The  third  curve  (.)  is  the  average  intra-si.^nal  standard 
deviation.  It  is  tvpicallv  a  nonotcnically  increasing 
function  of  signal  du'^ation  from  10C  ms  to  1^  seconds.  The 
fact  that  the  int ^'a-signal  standard  deviation  is  the  Sfriallest 
at  1'20  ms  indicates  that  the  major  factors  affecting  tne 
variance  ere  not  rapidly  fluctuating,  i.e.  the  period  is 
larger  than  I'id  ms.  The  intra-signal  standard  deviation  at 
100  ms  is  small,  four  to  eight  degrees,  except  for  file  1/5 
(?i,2;.  25-23)  which  was  recorded  with  a  low  SNR.  In  tnis  case. 
a  major  source  of  variance  was  noise  and  its  fluctuations 
were  rapid.  As  the  signal  duration  increases,  the  pnenomena 
causing  the  malority  of  the  variance  have  mere  effect  within 
individual  signals. 

The    fourth   category   of   data   on   the   graphs   is   the 
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probaoility  of  cttaining  a  line  of  tearing.  The  derivation  of 
this  number  has  previously  been  ex;)lained.  Its  usefulness  is 
a  proper  subject  for  operations  research;  it  does  not  provide 
the  engineer  with  any  useful  inf o rma tior  about  t ne  process  of 
HFDF  on  short  duration  signals.  However,  one  does  not  need  a 
specific  operational  context  to  know  that  hign  ?0B  is  good 
and  low  FOB  is  bad. 

One  can  summarize  tnis  section  by  stating:  (1'  Ai  is  an 
effective  measure  of  reliability  of  tne  data,  (2)  fie 
variance  is  high  for  short  duration  signals,  (3;  tne  average 
bearing  error  is  not  useful  without  calibration  data  and  (4) 
variance  improvement  with  time  averaging  corresponds  in  order 
of  magnitude  with  that  predicted  bv  fading  phenomena. 
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C.  LMAT 

The  previous  section  detailed  the  results  cf  accepting 
data  record  hearings  on  tne  hasis  of  the  Ai  ter.T;.  The 
calculated  variances  are  lar^^e  and  in  some  cases  the 
corresponding  PCB  is  small.  There  is  a  strong  need  to  lower 
the  variances  and  increase  the  POB  to  approximately  one.  To 
do  this  requires  a  more  complete  use  of  the  data;  tnis 
includes  the  k?'   and  Phase  terms. 

The  theoretical  equations  of  the  coaxial  spacec.  looris 
prove  the  eristence  of  polarization  independent  nulls  in  the 
ideal  case.  However,  the  construction  of  perfect  loops  and 
free  space  siting  of  these  loops  is  clearly  impossible. 
Pecognizing  that  the  spaced  loop  array  does  not  perform 
ideally,  one  attempts  to  identify  the  sources  of  fixed  error 
and  correct  them  with  a  calibration  data  set.  Otner  sources 
of  error  are  random  and  are  described  in  probabilistic  terms. 
It  is  usually  the  case  that  more  than  one  of  tne  major  errors 
is  random.  This  results  in  joint  probability  density 
functions  that  are  impossible  to  derive  analytically  and 
impossible  to  isolate  in  order  to  measure. 

In  an  attempt  to  maximize  the  use  of  information  available 
without  a  precise  knowledge  of  either  the  fixed  or  random 
errors,  it  was  decided  to  take  a  decision  theory  approach.  An 
assessment  of  a  bearing's  reliability  is  what  is  most  needed. 
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The   decision   to  te  made  is  tinaryJ  a  "bearing  is  reliable  or 
unreliable.  Once  the  tearin.^  reliability  is   deterniined,   tne 

statistical   procedure   developed   for   DFERR   can  oe  used  to 

determine  a  mean  bearing.  The   decision   to   determine   if  a 

bearing   is   reliable  will  be  based  on  the  li:;^elihood  that  it 
is  reliable. 

As  an  examDle,  assume  the  following  data  records: 


Record  #   A0 


A4 


bearing 


1 
2 
3 

4 


0  .113 
0.612 
0.514 
0.815 


0.214 
0.208 
0.421 
0.113 


340 
356 
266 
3^5 


The  true  bearing  is  icnown  to  be  337.   Assume   that   a   second 
data  set  is  recorded  on  an  unlocated  target  transmitter: 


Record  # 

A0 

A4 

Pear ing 

1 
2 

3 

4 

0.056 
0.822 
0.109 
'3.666 

0.298 
0.173 
0.202 
0.432 

105 
042 
093 
087 

Which  record  is  most  reliable?  It  is  known  in  data  set  ere 
that  the  first  record  is  tne  most  reliable  because  its 
"bearing  is  the  most  accurate.  Using  the  maximum  li^rcelinccd 
(ML)  criterion,  one  would  estimate  that  the  record  in  data 
set  two  that  most  closely  resembles  record  one  of  data  set 
one  is  the  most  reliable.  Bv  ML  c^'iterion  record  tnree  would 
be  selected  because  its  A0  and  A4  parameters  are  the  closest 
match  to  record  one  of  the  first  data   set.   In   this   scheme 
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there  is  no  filtering  on  preset  A0  or  A4:  limits.  If  tae 
lowest  level  of  k7>  and  A4  were  the  reliabilty  criterion, 
record  one  of  the  second  data  set  would  have  been  selected  as 
the  best  estimate.  The  advantage  of  the  ML  criterion  is  tnat 
it  is  a  "use  what  wor^^s"  technique.  The  deductive  approach  of 
analysis  of  the  system  and  errors  is  discarded  because  it  is 
too  complicated.  Instead,  the  inductive  tecnnique  of 
observing  and  classifying  provides  the  more  attractive 
approach  in  this  case. 

To  use  the  ML  criterion  with  the  SViRI  data  files  one  must 
expand  upon  the  ideas  presented  in  the  example.  Define  two 
matrices,  a  bearing  acceptable  and  a  bearing  unacceptable 
matrix.  Each  matrix  is  three  dimensional;  a  dimension  is 
allotted  to  A?',  one  to  Paase  and  the  third  to  A-i.  This  spans 
all  the  information  in  a  data  record.  A0  is  a  measure  of  tne 
horizontal  field  component;  Phase  is  a  measure  of 
polarization;  and,  A4  is  a  measure  of  system  inconsistency. 
These  three  measures  are  not  sufficient  to  completely  specify 
system  performance;  in  fact,  they  are  not  sufficient  to 
completely  specify  the  horizontal  field  component, 
polarization  or  syste*^  inconsistency.  Tney  are  what  is 
available,  "ach  dimension  is  divided  into  ten  increments.  The 
k?'  and  A4  dimension  values  are  between  2  and  1.-^  in 
increments  of  0.1  .  The  Phase  dimension  is  between  -1£2  and 
+  1S0  in  incT-ements  of   36   degrees.   The   structure   of   both 
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matrices    can    "be   sketched    as: 


PH^SE 


J — I — i — I — ' — I — I — J — « — I — 1 
O  0.5  ^'O 


PH^SE. 


~\80 


M 


0.5 


\80 


i.O 


The  matrices  are  defined  to  be  expressions  of  tne 
conditional  prohahility  mass.  They  are  created  by  using  a 
file  or  files  of  data  on  a  target  of  known  position.  Knovinf: 
the  true  bearing  permits  one  to  address  either  the  beariait; 
acceptable  or  bearing  unacceptable  matrix.  Further 
explanation  is  best  presented  using  an  example.  The  following 
records  are  available  from  a  data  set  with  a  true  bearing 
known  to  be  337  deg: 


Record  U      A 


AZ 

Phase 

A4 

Bearing 

(deg) 

(deg; 

0.911 

-43 

0.621 

286 

0.813 

-64 

0.518 

273 

0.517 

120 

0.231 

357 

0.342 

72 

0.165 

340 

2 

3 
4 


All    the   elements    of    the    bearing   acceptable   and 
unacceptable  matrices  are  initialized  at  zero.  The   criterion 

99 


for  placement  in  one  matrix  or  tne  other  is  the  value  of  the 
bearing.  A  window  aboiit  the  true  bearine-  is  defined.  If  tae 
window  is  defined  as  1^  degrees,  a  bearing  is  acceptable  if 
it  is  withi'^.  5  degrees  either  side  of  the  true  bearing.  "^acn 
record  in  a  file  is  examined.  Record  number  one  is 
unacceptable  because  the  beaming  value  is  not  in  the  332-3-±2 
deg  window.  Therefore,  its  AC,  Phase  and  A4  values  are  mapped 
into  increr^ents  along  the  respective  dimensions  of  tne 
unacceptable  matrix,  thereby  addressing:  one  of  the  elements 
of  the  matrix.  A  one  is  added  to  the  contents  of  tnis 
element.  Each  T'eccd  in  turn  is  examined  .and  a  one  is  added 
either  to  an  element  in  the  bearing  acceptable  or  the  bearing 
unacceptable  matrix.  (Record  number  four  is  an  example  of  a 
record  that  would  apply  to  the  bearing  acceptable  matrix.) 
After  all  records  have  been  processed,  each  element  in  the 
bearing  acceptable  matrix  is  divided  bv  the  number  of  records 
that  contributed  to  that  matrix.  3imila-^ly,  the  elements  of 
the  bearing  unacceptable  matrix  are  divided  by  the  number  of 
records  that  contributed  to  it.  This  produces  an  expression 
of  conditional  probability.  If  a  bearing  is  acceptable,  the 
probability  of  it  having  a  particular  A^,  Phase  and  A-i  value 
is  equal  to  the  value  in  the  acceptable  matrix  addressed  by 
the  given  A0 ,  Phase  and  A4:  values.  Likewise,  if  a  bearing  is 
unacceptable,  the  probabilities  of  its  Ad,  Phase  and  A4 
values  can  be   read   fT'om   the   element   addressed   by   those 
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values  . 

Tne  acceptable  and  unacceptable  natrices  are  constructed 
fron  data  on  a  located  target.  If  one  Drocesses  a  file  of 
data  en  an  unrcnown  tar^^et,  the  matrices  are  used  as  follov/s. 
Using  the  A?,  Phase  and  A4  values  to  address  both  of  the 
matrices,  the  element  values  of  the  matrices  are  compared. 
Supp'^se  the  value  from  the  acceptable  matrix  is  '^.211,  and 
from  the  unacceptable  mat-ix,  it  is  2.^9?.  The  ML  criterion 
makes  the  decisio'^  fo^  the  matrix  with  the  aignest 
probability  (maximum  li 'rcel  ihood  )  .  ^.'Sll  is  the  highest  value; 
therefore,  the  decision  is  that  the  bearing  is  acceptable. 

The  requirement  to  address  two  matrices  and  compare  tae 
returned  values  can  be  eliminated  by  forming  a  single 
likelihood  ratio  matrix.  This  matrix  is  constructed  by 
dividing  each  element  of  the  bearing  acceptable  matrix  by  the 
corresponding  element  of  the  bearing  unacceptable  matrix.  The 
decision  can  be  made  on  a  single  element  of  tne  likelihood 
ratio  matrix  (also  simply  called  likelinood  matrix  or  denoted 
as  [L]  )  .  In  the  example  above,  the  ratio  of  the  two  values  is 
2.175.  The  decision  is  that  a  bearin.?  is  acceptable  if  the 
addressed  element  of  the  likelihood  matrix  is  greater  than  or 
equal  to  -^re.  If  less  than  one,  the  bearing  is  unacceptable. 

The  above  ideas  ca-^  be  set  into  mathematical  notation  as 
follows,  i'^at^ices  ave  denoted  witi  brackets.  The  elements  of 
a  matrix  are  indexed  with  i ,  .i , -^  in  the  general  case. 
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[A.]  t  bearing  acceptable  matrix 

f^"^]  :  "bearing  unacceptable  rnat^'ix 

[A    (i.^.k)    ]    =    Ta    (  A0,  Phase,  A4:)1 

fn    (i,j/-c)    1    =    (U    ^A0, Phase  ,A4H 

Na   =   number    of    records    tabulated    in    lA] 

Mu    =   nurnbe^    of    records    tabulated    in    [U.i 

ila    +    Mu    =    {l'd'dd(3)    X    (number    of    files    used 

The    DTobability  mass   matrices   are: 


[Pal    =    (1/Na) lA] 
(?ul    =    (1/Nu) fU] 


The  maximum  lilcelihood   fML)   criterion   can   be   used   to 
decide  the  acceptability  of  a  beaming. 
Let  the  record  be  f A3=^^,?hase=A.A4='<.'Bo } ,  then. 


DZ. 
> 


"^cA"^^    0\ 


TT 


^c/?.n 


where,  the  decisions  Dl  and  D2  are. 

D*^. :  Bearins;  is  acceptable 
Dl :  Bearing  is  unacceptable 

If  A  >  U.  the  decision  is  J)?'    if  A  <   U,  the  decision  is  Dl. 

This  compa'^isor  process  car  be  simplified  and  at  tne  same 
time  be  made  more  flerible  by  defining  a  liicelinood  ratio 
matrix  [L"!  : 

fLl  =  fPal  ©  [^u] 


1?2 


where  the  symbol  (©)  means  to  divide  each  element  of  the 
matrix  [Pa]  only  by  its  i ,  j ,  k:  counterpart  in  '?ul.  (Define 
T/'l    =    infinity.) 

The  ML  criterion  can  be  rewritten  as: 


[LI 


Da 

> 
< 


[1] 


[1]  =  13x12x1^  matrix  of  1 's 


D2 
> 


The      above    idea    can   be    extended    to    the    I'^inimum   Probability    of 
Error        ("-^P^J        and        the        Saves  Cost  (BZ)  criteria 

s  traiiSihtf  orv/ardly : 


D2 
fLl   ^ 

Dl 


Ll 


D2 
> 
< 

1)1 


P(3  acceptable} 
PiB  unacceptable} 


Cm  • 

^   12.  • 

C  22.  : 


(C^,  -  C,,  )     PiB  ace} 
^^  \-L    ~    ^  12.^  '^-^   unacc} 


cost  of  deciding  acceptable  wnen  unacceptable 

cost  of  deciding  unacceptable  wnen  acceptable 

cost  of  decidinf^  unacceptable  when  uiacceptable 

cost  of  decidin.^  acceptable  when  acceptable 


The   terms   on   the   ri^ht   side   of  the  decision  symcols  are 


/   matrices.  The  ratio 


P{3acceptable}   to   P(& unacceptable} 


mi:st  be  determined  separately  for  each  element  of  the  matriy. 
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The  sane  will  probatly  "be  true  for  the  Ba/esian  costs.  Tne 
MPE  and  5C  criteria  are  nore  sophisticated  and  will  protatiy 
yield  better  results,  but  data  limitation  prevented  further 
investigation.  Very  large  amounts  of  data  would  be  needed  to 
determine  probability  ratios  for  each  element  of  a  10xl<5xl3 
matrix.  It  would  also  require  a  operational  input  for  a 
realistic  assessment  of  3avesian  costs. 

Using  the  above  concepts,  the  FORTRAN  program  LMAT  was 
written  to  produce  the  li'Celinood  ^atio  matrices  [Lj  f-^om 
individual  files  and  from  groups  of  files.  Testing  of  tne 
method  centered  on  the  15  "^Hz  files  because  they  are  the  most 
numerous  for  a  single  frequency  (three  files).  An  L  matrix 
was  created  from  file  3  for  two  separate  definitions  '^f 
"bearing  acceptable".  The  first,  denoted  Li'3,9'?J.  considered 
the  bea^'ing  acceptable  window  to  be  ninety  degrees  wide.  The 
second,  L^Z.ZTi],  used  a  twenty  deg^^ee  window  to  construct  trie 
matrix . 

Using  the  L  matrices,  several  data  files  were  processed 
employing  the  program  LFILE  (L  matrix  modified  FILji^j.  Tnis 
program  reads  the  A0 ,  Phase  and  A4  term  of  each  record  and 
uses  them  to  address  the  L  matrix.  If  the  element  addressed 
is  greater  than  one,  the  bearing  in  tnat  record  is  considered 
acceptable,  and  the  record  is  written  into  a  new  data  file 
with  a  new  A4  term  equal  to  ?.?1.  Onlv  the  A4  term  is 
changed.  The  new  data  file  is  subsequently  processed  ty   DF5RR 
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which  will  always  treat  an  acceptable  bearing  {ki=d .dl)  as 
valid  for  statistical  processing.  If  the  L  matrix  element 
addressed  is  less  that  one,  the  bearing  is  unacceptable  and 
assi.-?ned  the  A4  value  10.!^  in  the  new  aata  file.  DFiRH  will 
reject  any   bearing  with  this  lar^e  A4  value. 

The  file  of  interest  is  the  1990  15  MHz  file.  ?i.-^ure  35  is 
the  result  of  simple  DF?]?.R  processing  on  this  file;  figure  56 
is  the  hlsto^^am  of  the  bea'^ing  e'^ror.  Several  examples  of  L 
matrix  processing  nn  this  file  a^e  .-graphed  in  fi=:ures  37,  c=j 
and  4:1.  The  explanation  of  these  e-r'aphs  is  identical  to  the 
explanation  ^iven  in  the  previous  section.  It  will  be  noted 
, however,  that  the  graphs  d.-re  labelled  with  the  L  matrix 
notation.  The  A4!^AX  term  is  not  applicable.  Each  graph  has  an 
accompaning  histogram  o^  the  b*=»aring  er^^o-'.  Figure  Z'^  is  tie 
file  processed  by  an  L  matrix  made  from  tne  file  itself. 
Figure  39  is  the  file  processed  bv  an  L  matrix  created  from  a 
1979  file  at  the  same  frequer'^y.  Figu^'e  41  is  the  file 
processed  bv  a^  L  matT-ix  created  f^om  fou^  1979  a.vV  files  at 
different  frequencies.  Mot  unexpectedly,  the  best  performance 
is  by  the  L  matrix  formed  from  the  file  itself  (Fig.  37).  But 
it  is  interesting  to  note  that  the  L  matrix  created  from  the 
single  file  of  identical  freauency  (Fig.  39)  performed  better 
than  the  matrix  made  up  from  the  four  files  (Fig.  41).  This 
tends  to  confirm  frequency  sensitivity  in  the  L  matrix.  It 
may  prove  valuable  to  add  a  frequency  dimension,  making  the  L 
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matrix  four  dimensional. 

A  considerable  amount  of  processing  was  devoted  to 
determine  an  optimum  vindow  for  the  likelihood  matrix.  No  one 
window  width  outperformed  other  widths  in  all  categories.  The 
central  difficulty  is  that  one  thousand  elements  in  the  L 
matrix  must  be  determined.  Accuracv  in  terms  of  i^L  for  each 
element  is  a  limiting  process;  as  the  number  of  bearin=^s  used 
to  deter^iine  each  element  apprcacnes  infinity,  tne  true  i^^L 
ratio  is  determined.  To  approximate  tne  infinite  sample  ML 
ratio  within  1'^  "oercent  would  require  very  rougnly  120 
records  per  element.  But  some  elements  are  rarely  adaressedj 
therefore.  large  amounts  of  data  are  needea  to  fill  tne  L 
matrix.  The  amount  of  data  can  be  estimated  in  a  very  ^ough 
way.  Usin^  four  ^iles  and  a  ninety  decree  window,  it  was 
obseT'ved  that  ^bout  40  percent  of  the  matrix  elements  a.re 
addressed  zero  or  one  times.  If  the  requirement  is  tnat  the 
probability  of  this  happening  should  be  less  than  0.201.  tne 
expression  (0.4)  =(0.001)  determines  the  number  (n)  of  sets 
of  four  files  ^equi^^ed.  In  tnis  case  about  thirty  files  would 
be  required.  This  corresponds  to  300,000  records  or  100 
minutes  of  data.  This  is  an  easily  achieved  number  at  the 
antenna  system  devel  ODme'i  t  site. 

The  most  overall  successful  runs  were  made  with  a  ninety 
degree  window.  In  table  VI  are  the  results  of  processing  tne 
1960  data  files  with  a  likelihood  matrix  constructed  from  the 
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1979  data.  The  first  four  1979  files  were  use'i  ;  the  fiftn,  a 
second  15  MH^  file,  was  not  included  to  avoid  a  possible 
dominant  15  '^Hz  bias.  The  data  in  table  VI  is  the  variance 
and  P05  from  sinple  D?EF?  processing  ( A4:MAX-=0 .2  )  and  fT'on  L 
matrix  processing  (L  ( l-4:,90 )  ■  .  For  the  193-^  5  MHz  file  the 
variance  is  high  in  either  case.  Starting  at  30^  ms  signal 
duration,  the  L  r^atrix  processing  has  a  lower  variance,  a.id 
as  the  signal  duration  increases,  L  matrix  processing  is 
better  and  tetter  conpared  to  DFEHR  processing.  Pernaps  mucn 
"lore  significant  is  that  the  L  processing  has  a  nuch  nigher 
PCE.  .Vith  the  19EZ  1?  MHz  file.  L  matrix  Drocessing  is 
superior  at  everv  signal  duration.  This  is  also  true  i or  the 
15  MH7  file.  This  is  a  significant  result.  L  matrix 
processing  is  generally  suT^ericr  tc  the  A4  filter  process  cf 
DFZRR.  In  sone  cases  the  variance  is  halved  and  tne  PC3  is 
more  than  doubled.  HoveveT',  performance  against  iLC  is  poorer 
as  can  be  seen  in  the  table.  For  tne  shortest  signal 
durations.  the  L  matrix  variance  is  higher  tnan  DFEPH 
processing,  but  tne  two  approach  each  other  past  500  ms.  Tne 
FOB  is  slightly  higher  for  the  L  matrix  process. 

The  results  in  table  VI  are  very  encoura=?ing.  A  matrix 
created  from  1979  data  significantly  outperformea  DFERR 
pT'ocessing  on  the  198?  data.  Tne  poorer  results  for  fJLC  are 
not  surprising.  'a-.W  bears  337  deg  from  San  Antonio,  a^d  K^C 
bears  ?9?   deg.  The  antenna  array  which   is   mounted   atcD   a 
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hcri  zcn  tall""'  girded  "'buildir.g  suffers  patterr  aistortion  due 
to  antenna  eler^ent  co-'jplinf:  with  the  "building.  Tne  distortion 
to  the  pattern  is  dependent  both  on  the  azimuth  and  elevation 
angle  at  which  the  signal  arrives.  In  fact,  tne  KLC  results 
are  encouraging.  Tne  fact  that  L  matrix  processin<?  and  DFERR 
processing  results  are  so  close  suggests  tnat  azinuthal 
dependence  is  not  an  overiding  factor.  It  rray  be  possible  to 
achieve  acceptable  i^esults   usin>?   L   'matrices   created   fro.^ 


azir^uth  sector  data.  I 


^i 


irst  guess  is  tnat  a  bearing  is  in 


the  sector  21?  to  27''^,  post  processing  could  use  an  L  ratr'ix 
specifically  ceated  for  that  6<?  degree  sector.  In  all  six  L 
matrices  would  be  reqvired  to  cove"  the  full  aziruta. 

The  optimism  expressed  above  must  be  tempered  by  the  fact 
that  the  technique  and  the  results  were  derived  with  an 
inductive  approach.  A  theoretical  basis  with  wnich  tn.e  L 
matrix  success  could  have  been  predicted  was  not  derived.  L 
matrix  success  was  observed,  not  predicted.  Confidence  that 
the  technique  is  functional  in  the  general  case  will  require 
observation  in  situations  whe^'e  f-^equency.  azimutn.  3NK, 
vertical  angle  of  aT'rival,  polarization,  shipboard  siting  and 
propagation  modes  a^'e  va^^ied. 
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TABLE    71 

Comparison    of    DFEP.P   and    L   '"Matrix   Processin.^ 
?or      all      DPER?      iDT'ocess  inf: ,      A4:MAX  =  0 .2 .      For      all      L      matrix 
processing,    the   L   rat^ix   is    1(1-4. Q'J). 


5if?nal    Duration    in   millis  econds  : 

1?0      2'^0      Si?-?      4:V?0      o?0      6c'0      702      500      900       1000 


W'iV  5 

MHz 

20  :40 

2/6/ 

30 

3TD 

DfFP? 

57 

59 

5Q 

59 

59 

63 

61 

o6 

62 

55 

LMAT 

62 

62 

5i 

49 

46 

40 

40 

36 

36 

o5 

?0B 

DFE^?. 

.11 

.33 

.52 

.  65 

.77 

.83 

Q  /I 

.91 

.95 

.93 

LMAT 

.04: 

.93 

.99 

1 

1 

1 

1 

1 

X 

1 

1 

WyV  10  MHz 

7:40 

2/5/ 

80 

STD 

DFERR 

30 

29 

25 

21 

13 

17 

13 

16 

10 

11 

LMAT 

29 

23 

16 

10 

11 

c 

7 

9 

5 

6 

FOB 

DFERP 

.73 

.90 

.95 

.97 

.99 

.99 

1 

1 

1 

1 

LMAT 

.90 

.97 

.99 

1 

1 

1 

1 

1 

1 

1 

>/;^V  15  MHz 

9:00 

2/5/80 

STD 

DFERR 

16 

16 

14 

12 

13 

11 

9 

9 

9 

9 

LMAT 

13 

13 

13 

10 

9 

10 

9 

10 

10 

9 

?05 

DFERR 

.P3 

.86 

.87 

.88 

.89 

.89 

.9? 

.90 

.91 

LMAT 

.ai 

.86 

.88 

.90 

.91 

.91 

.92 

.93 

.94 

KLC  8.^ 

':e6    MHz  8  :: 

15  2/ 

7/30 

STD 

DFERR 

1^ 

19 

16 

16 

14 

14 

11 

15 

10 

10 

LMAT 

23 

20 

22 

20 

18 

17 

12 

11 

li 

9 

FOB 

DFERR 

.22 

.^d 

.43 

.46 

.47 

.43 

.50 

.50 

.54 

LMAT 

.*25 

.42 

.46 

.43 

.48 

.50 

.t^Z 

.52 

.56 
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D.  AMBIGUITY  Ri-SOLUTIOM 

The  his  tOf^T^ans  of  the  ave^a^e  bearing  e^'^or  demonstrate 
one  of  the  difficulties  vlth  trie  3WFI  spaced  loop  antenna 
syst^n  in  t'^yinp  to  DF  sho^t  duration  signals.  It  can  "be  seen 
cy  loo'-cin^  at  the  histoprans  that  the^e  are  "bearing  error 
numbers  that  cluster  in  relatively  lar^e  values  at  points  of 
1S0  and  90  degree  ambi^iuity.  This  corresponds  to  tie  antenna 
system  deciding  on  the  vron^  null  in  tne  antenna  array 
pattern.  The  severity  with  wnich  this  can  affect  the 
calculated  variance  is  graphed  in  fi2;ure  43.  This  graph, 
plotted  f^op  an  equation  developed  in  a:)pendix  3. 
demons  t ''a  tes  that  a  small  number  of  ambiguities  have  a  ve'-y 
significant  effect  on  the  variance.  The  graph  assumes  1020 
bearings  constitute  the  total  sample.  All  the  ambiguities  a^e 
at  183  degrees.  The  abscissa  is  the  number  of  13f;  degree 
ambiguities.  The  ordinate  is  the  ratio  of  the  nev  variance  to 
the  old  variance  which  is  set  at  twenty  degrees  squared.  If, 
for  example,  there  are  ten  190  degree  ambiguities  in  tne 
original  samDle  and  they  are  removed,  the  new  standard 
deviation  will  be  less  than  ten  degrees  squared.  In  the  case 
of  many  small  sets  representing  short  duration  signals, 
evaluated  separately.  the  reduction  in  variance  will  be 
larger  than  that  p^^edicted  in  figu^'e  43. 

To  examine  the  effect  of  removing  the  ambiguities  from  the 
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SVRI  data  files,  a  FORTRAN  prograni  A^^BIa  was  written  to  set 
the  A4:  value  to  Id.d  for  every  record  wnich  is  within  a 
narrow  window  about  the  9?  and  1S0  degree  anoiguities. 
Setting?  the  A4  to  such  a  large  value  perr^itted  the  use  of  tne 
DFERR  program  which  automatically  rejects  records  witn  such  a 
high  A4  ter'".  This  test  was  run  on  several  data  files  and  the 
results  are  displayed  in  figures  44  and  45.  It  can  be  seen 
that  snarp  decreases  in  variance  occur,  serving  to  Tia.-ce  tne 
bearing  estimate  more  reliable. 

The  resolution  of  ambiguities  with  narrow  apert^i^e  syste^^s 
is  a  difficult  p^'otlem.  The  alro^'ithms  are  especially 
sensitive  to  vertical  angle  of  a^ri'"'al,  low  SMR  and  m.ultinode 
propagation.  ""ime  averaging,  if  the  signal  is  sufficiently 
long,  can  overcome  multimode  propagation,  but  vertical  angle 
of  arrival  is  very  dependent  on  array  geometry,  and  S.MR  is 
often  totally  uncontrollable.  The  spaced  loop  array  amDiguity 
problem  may  be  due  to  equipment  errors,  but  it  is  more  likely 
that  the  physical  limitation  of  the  aperture  is  the  primary 
difficulty.  The  aperture  is  insufficient  to  sense  wavefroat 
distortions  or  to  sharply  isolate  pattern  nulls.  In  tne  case 
of  burst  communications,  time  avera.^ing  is  impossible.  The 
narrow  aperture  array  must  be  assisted  by  wide  aperture 
elements  to  work  against  short  duration  signals.  If  a  spaced 
loop  array  is  "^ounted  aboard  a  naval  vessel.  it  may  be 
possible   to   add   a  simple  loop  el  erne':  t  at  the  bov;  and  stern 
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and.  two  amidships.  This  sij^ple  cross  shaped  interferometer 
array  would  be  medium  aperture.  By  itself  it  would  not 
provide  reliable  H?DF,  but  it  could  improve  the  reliability 
of  tne  narrow  aperture  system  by  sensin^^  wavefront 
distortions  and  providina:  deep,  well  defined  nulls  to  resolve 
ambi^ui  ties . 
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Ambiguity  Suppression  Curve 
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Figure  44 
Ambiguity  Suppression  for  WWV  5  MHz  2/80 
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Figure  45 
Ambiguity  Suppression  For  WWV  15  MHz  2/80 
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VI.  CONCLUSIONS  AMD  P  ilCOi^MINDATI  QMS 


A.  CONCLUSIONS 

The  narrow  aDerture  antenna  is  Dhvsicall.y  incapatle  cf 
spatially  sensin.^  phase  or  amplitude  wavefront  distortions 
that  have  spatial  periods  of  manv  wavelen,^tns .  This  inherent 
disadvanta<^e  ^^ust  be  overcome  ty  tine  averas:in<2:.  Under  t;ie 
assumption  that  wavefront  distortions  fluctuate  ahcut  tne 
undistorted  wavefront  (the  mean  wavefront),  the  minimum  time 
necessa'^y  to  average  out  the  distortions  is  equal  to  the 
lon,sest  period  of  the  primary  pnenomena  causin,^  the 
distortion.  The  most  severe  distortions  occur  with  multimode 
interference  when  two  or  more  rays  are  comparaole  in 
amplitude.  In  this  case  the  time  required  to  average  out  the 
fluctuations  is  the  fadin,^  period  of  the  major  spectral 
component  of  the  fading  power  spectrum,  ro^i^hly  on  the  order 
of  ten  seconds  for  an  HF  signal .  This  means  that  a  narrow 
aperture  antenna  cannot  ^eliatly  determine  the  direction  of 
arrival  of  a  hurst  transmission  in  the  presence  of  severe 
multimode  interference.  As  the  seve'^ity  of  distortion  is 
reduced,  the  reliability  of  D?  is  imprrved.  In  tne  case  of 
the  SV/kI  spaced  loop   fIFD?   system,   the   range   of   standard 
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deviation   for   the   200   ms   signals  studied  here  is  from  15 

degrees  for  signals  sysnthesized  from  data  in  file   4:   to   59 

degrees   for   signals   synthesized   from  data  in  file  1'.?.  Tp.e 

average  standard  deviatif^n  over  the  nine  data  files   is   29  1 
degrees.   These  data   files   a ^e   Relieved   to   include  ooth 

multim'^de  and  single  mode  signals  and  si^-'nals  with   ni^h   and 
low  SNF's. 

2.  Performance  Specifications 

The  performance  of  a  nar'row  aperture  ESD?  system 
against  short  duration  signals  is  very  dependent  on 
ionospheric  propagation.  Therefore.  the  performance  of  t.ie 
system  cannot  he  specified  separately  fr-^m  specifications  on 
the  state  of  the  ionosphere.  To  state  that  a  svstem  must 
perform  T</ith  a  given  standard  deviation  of  Dealing  error 
during  quiet  ionospheric  conditions  with  single  mode 
propagation  is  mucn  different  from  requiring  the  sar^e 
standard  deviation  of  bearing  error  during  disturbed 
ionospheric  conditions  or  with  multimode  propagation. 


5  •  Z^^iHC  §11^,  Sea  ri  ng  ^Piiah i  1  _i  t  v 


There  are  several  sources 


fading   in   ionospheric 


propagation;  however,  fading  can  be  considered  a  good  measure 
of  the  distortion  of  local  constant  phase  and  constant 
amplitude  wavefronts  .  If  fading  is  severe,  on  tne  order  of  12 
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d"b ,  one  can   be   confident   that  sL^nificant   distortion   is 

occurring:   to   the  vavef^ont  and  t!iat  narrow  aperture  derived 

DF  bearings  on   sho^t   duration  signals   will   be   generally 
unreliable . 

Tne  concept  of  using  a  iikelinood  .matrix  develoned 
from  given,  reliable  data  to  ral-:e  a  binary  (acceptable  or  not 
acceptable)  decision  on  a  random  data  set  has  not  been  fully 
explored.  However,  the  results  from  the  limited  processing 
are  very  encouraging.  It  should  prove  possible  to  develop 
appropriately  sized  and  numbered  L  matrices  capable  of 
significantly  outuerf orming  algorithms  based  only  on 
filtering  bv  laarameter  limits. 

5.  AlD^iSliiAZ  Elimination 

The  filtering  to  eliminate  Q')  and  1S0  degree 
ambiguities  demonstrated  tnat  a  significant  reduction  in 
bearing  va^ia'^ce  can  be  accomplished  if  the  source  of  the 
ambiguities  can  be  corrected.  It  is  doubtful,  however,  that 
the  correction  can  be  accomplished  without  the  introduction 
of  supportive  medium  aperture  interferometer  elements. 
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B.  RECOMMENDATIONS 

!•  Use  of  All  Ionospheric  Data 

HFDF  with  narrow  aperture  antennas  on  short  duration 
5i,2;nals  rust  maTimize  the  use  of  as  much  real  tire 
information  as  possible  on  the  ionosphe'^e.  Evaluation  of  the 
vertical  an^^le  of  arrival  and  polarization  is  needed  to  help 
determine  if  a  ground  wave  or  a  sk:vwave  tu^st  transmission  is 
heing  received.  The  vertical  angle  of  arrival  will  further 
permit  an  estimate  of  range.  The  range  and  vertical  angle  of 
arrival  information  should  he  used  with  a  propagation 
prediction  program  to  make  an  overall  evaluation  of  tne 
reliatility  of  a  calculated  bearin,?  based  on  predicted 
ionospheric  induced  tea^'ing  variance.  To  improve  tne 
propagation  prediction  p^o^ram,  ■^eal  time  measurements  of 
ionospheric  parameters  are  needed.  This  could  be  accomplished 
by  updated  inputs  of  geophysical  data,  shipboard  ionospheric 
sounders,  solar  observations,  data  derived  from  satellite 
beacons  and  the  use  of  geographically  fixed  transmitters  as 
beacons.  The  provision  of  this  information  will  probably  not 
permit  the  calculation  of  a  more  accurate  bearing,  but  it 
will  help  to  determine  the  reliability  of  the  bearing  (the 
variance  of  the  bearing). 
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2-  n^e  2f  All  Signal  Parametej:  Data 

Various  receiver  parameters.  especially  the  AGC 
volta<^e.  should  be  monitored  to  estimate  the  fading  of  a 
sis;nal.  A  lar^^e  variation  of  the  AGC  voltage  could  sis-nify 
fadin=^  and  somewhat  reduce  the  confidence  in  the  calculated 
tearing.  A  small  variation  in  the  AGC  voltage  would  be 
inconc lusi ve . 

?•  P?^5i.i§l  I§l£ormance  Specifications 

The  specifications  for  the  reqiiired  performance  of  a 
narrow  aperture  antenna  system  must  incliide  the  ionospaeric 
conditions  under  which  performance  is  to  be  measured.  An 
example  specification^  is  that  a  system  must  have  ^^n  ave^a^e 
bearing  error  of  five  lies^ees  and  a  standard  deviation  that 
encompasses  at  least  67  percent  of  the  data  when  tne  signal 
is  received  in  the  20  meter  band  via  a  predominan tl v  one  nop 
propagation  path.  The  signal  power  of  other  pT'oparation  modes 
should  be  30  db  bel'^w  the  primary  mode,  and  tne  SNR  against 
backgrni:nd  noise  should  be  at  least  12  db  .  To  conduct  such 
measurements  aboard  a  ship  may  require  that  the  ship  be 
positioned  close  to  a  wide  apertu^^e  array  that  would  be  able 
to  resolve  the  local  node  structure  of  a  tcnown  target 
transmitter's  signal. 
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4.  Likelinood  Ratio  Mc'.trix  Techniaue 


Though  the  results  of  usins  the  pro ha oali s t ic 
likelihooi  natrix  vrere  were  not  totally  conclusive,  it  is 
reco^nmended  that  this  approach  ce  further  inves  t  i^*a  ted  .  The 
stroTie^  3Dpeal  of  the  likelihood  matrix  approach  is  the 
maxi'^um  use  of  information  in  ar  i"^prec  isely  known  and 
inprecisel;''  knowable  environment.  Furtner  invest  ir:a  ti  on 
should  include  the  following  areas.  Optimize  the  dimensions 
of  the  likelihood  matrix,  perhaps  eliminating  tne  phase 
dimension.  Define  the  likelihood  matrix  in  terms  of  tne 
minimum  probability  of  error  criterion  insteaa  of  the  maximum 
likelihood  criterion.  Define  the  likelihood  matrix  in  terms 
of  Bayesian  costs  and  compare  tne  resrlts  with  those  of 
maximum  likelihood  and  minimum  probability  of  error. 
Determine  tne  frequency,  3NR  and  fading  sensitivies  of  the 
likelihood  matrix.  Of  greatest  value  and  greatest  difficulty 
would  be  a  general  theoretical  development  of  the  liKelinocd 
matrix  in  te'^ms  of  tne  spaced  looo  system  parameters.  It  may 
prove  to  be  the  case  that  tne  likelinood  matrix  could 
concisely  store  calibration  data.  This  could  be  calibration 
in  azimuth,  elevation,  polarization,  frequency  and  5Nr.,  all 
in  a  finite  number  of  likelihood  matrices.  Using  the 
measurements  of  frequency  and  SNP  a^A  the  estimates  of 
azimuth,  elevation  and  polarization,  a  stored  likelihcrd 
matrix  would  postijrocess  s^^stem  data  to   refine   tne   d earing 
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es  tima  te . 

5«  Apti£]ait,jr  Diagnostic  Algo;[;i  thm 

There  is  a  need  to  develop  a  dia.^nostic  alATcritnm  that 
will  use  a  '-cnown.  fixed  transmitter  to  tabulate  ambiguity 
error's  as  a  check:  o-^  system  sens  iti  vi  t:'  and  of  the  phase 
sensing  elements,  par ticularily  the  phasenete^. 

The  spaced  Iood  array  does  not  provide  reliable  data 
for  various  important  tactical  situations.  The  variance 
associated  with  a  burst  signal  bearing  is  too  large  for  most 
fix  and  targeting  algorithms.  To  imT)rove  performance,  the 
feasibility  of  adding  a  si-^Dle  medium  aDerture  interferometer 
should  be  investigated. 
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APF3NDIX  A. 


AViRAGil  AND  STANDARD  DEVIATION  CALULATiCNS  ON  BSAEL.GS 

The  con^piiter  pT'Oprarns  used  in  tne  analysis  portion  of  tnis 
report,  corputes  avera£;e  and  stania^d  ieviatio:".  of  cearin^-s. 
These  statistical  caloulat  io'^s  Are  s  t  ^aip-htf  o^va^d  ,  cut  not 
in  the  forn  (generally  reco^rnizei.  In  this  report.  teai^ir^TS 
have  been  treated  as  integers  in  the  set  3  to  3by.  When 
calculating  averages  of  tearings  that  overlap  the  ^^-559 
toundary.  it  nust  te  taken  into  account  that  359  differs  fror. 
0  by  only  one  degree.  For  this  reason  the  lolloving  fornulas 
were  ieveloDed. 

Assu'^e  that  n  bearings  are  available  for  averaging: 
Bl,62 3i Bn 

Nornally  the  average  would  be  calculated  by: 

B  =  (1  'n)  X)  3i 
Fewrite  the  bearings  in  terms  of  a  reference  bearing.  B  ref 
and  a  difference  Ai. 

Bi  =  B  ref  ^    Ai 


Then  . 


E  =  (1/n)   Z  'B  ref  ^  Ai)  =  B  -ef  +  (l/'n  'X!  Ai=  B  ref  -A 
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Average  bearing  =  Reference  3earin.=!:  +  Avera.ge  Difference 

It  is  nest  convenient  to  let  B  ref  =  '3j    then, 
B  =   A 

However,  the  r?reate5t  difference  allowe':^.  in  a  closed  set  ?  to 
359  is  If?;  therefore,  pax  =  180,  and  the  conditional 
definition  is: 

Ai  =  Bi   if   Bi  ^  ISe 

Ai  =  Pi  -  36  3   if   Bi  >  13  0 

Now  the  avera^^p  f'^'^'^ula  can  "be  written: 

B  =  (l-'n)E  A  i 

subject  to  the  conditional  definition  of  ,  wnic.d  is  easily 
inplenented  in  FOPT^AM  . 

Por  the  standard  deviation,  a  typical  formila  is: 


^B  =   'l/n)E  3i 


(B) 


Af^ain  let  : 


Ai  =  3i    ^  <  Bi  <  160 

Ai  =  £i  -  350     190  <  3i  <-  35,- 


Then,  in  a  lirce  ranne^. 


<^C=      (1/r.^EAi   -((l/n)EAi/^ 
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APPENDIX  3 


AMBIGUITY  SUPPRESSION  VARIANCE  EQUATION 

In  the  analysis  portion  of  this  report,  the  results  of 
supDressi  n<?  9.-5  and  16(^  de.^ree  anbie^ui  ties  from  a  "tearing  data 
set  were  discussed.  The  equation  used  to  descrite  tne 
significance  of  anhi^nity  suppression  is  derived  in  tnis 
appendix  . 

Assume  that  a  his  to.;; ram  of  bearing  errors  closely 
approximates  the  Trohabil t iy  density  function  p^b/.  where 
p(b)  is  sketched  as: 


p/z 

4 


4 


i 


P/a 

4 


^b 


-180 


-0(0 


^O 


180 


The  '^elta  functions  rep^^esent  that  po^tinn  of  tne  histogram 
due  to  ambi^^ui ties .  The  function  ^,{t)  is  j^ene^al;  tne  only 
restrictior  associated  with  it  a'^e  that  the  o-"igin  of  the  b 
axis  coincide  with  the  mean  of  £^(b^  and  that  b'b)  be  wide 
compared  with  the  distribution  of  the  ambi^'ui  ties . 

One  writes  ^(bi  as: 
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p(b)    =   ^'b)    +    ^?/2)S(b^l£0)    +   LS(b+9?)    ■'■   ?.  S(b-9v^ 

Usin^       /pCb)    db   =   1 

and  Afbldt   +    P   +    L   +   P^l 

trien  A(b)    db   =    l-{P+L    +   R)  (2) 

The   variance    of   "d(d^    is: 


VARfB]    =/*b^p(b)    db      -        /b 


TD  (  b  )    db 


/b   TD(b)    db   =   /b   ^(b)    db   -^    (P)(18'^)^    +    (L  +  P. '(9^;^ 


/ 


b    ID  (  b  ^    d  fc   =    d   +    (  L  }  (  -9^  )    +    i  P. )  (  9^^ 


-L;  (90) 


=/.% 


VARlBl     =/b    ^(b)    dt    +    (?  )  ( 180  )+ (  L+R  )  (93  '-(  ~-L  ;     (90)       (5) 


If    the    arr'bi^ui  t  ies    ai^e   eliminated,    a    nev   pdf    tT''vb)    results 


— 1 — 

-HO 


^W 


*io 


g'Cb)    is  a  scaled  version  of  i°:(c');  therefore, 
K  ^'(b)  =  g(b)   where  K   <  1 

)  db  =  1   ==>  / .?(  b;  dh  =  K 

copihinins   tnis    with    (2)    yields: 

K=1-(P+L+R)  (4) 


/-(V 


Tne   variance    of    the   new   T^df    is: 
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VAPfS']  =  /b  s'{b)  db   -   /b  ^'(b^  a 

ft'^'i^^    db  =  fl'I^  Tb^^'^b)  db 

/b  g'fb)  db  =  '3   (by  assumption,  centered  at  the  ori.c^in 
VAE  [?n  =  (1/K)  A^frlb)  db     (5^ 

Substitute  (5)  into  (3)  3nd  solve  fc  the  new  variance; 


VARI'E]  =  K  VAP.rs'l  ^  (?)(190)   ^  (L^R)(9Z'   -  (P-L)i9i'; 


VAP'P'l  - 


VAR  [B ]  ^  (  ( R-L)  -  ( L+R )  -  4?  M  90 ) 


I 


1  -  (?  +  L  +  R) 
Equation  (5'  is  the  expression  that  cDprcxinates  tne  new 
variance  of  a  general  districution  wnen  the  ambiguities  are 
eliminated.  Taere  is  one  condition  that  -nust  be  observed. 
Because  variance  is  always  a  positive  quantity,  the  numerator 
of  (6'|  must  be  a  positive  quantity.  (The  denominator  is 
quaranteed  to  be  positive.)  Therefore, 

2.  i 

VARfB]    +    ((R-L)    -      (L+R)      -      4P)(90)       >      0      wnicn      Con       be 


re-expressed    as , 
^g  >    9'd  14? 


+    (L^R^    -    LR 


,-.'] 


(7) 
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APPENDIX  G 


PROP  HZ T 


The  data  used  to  coripile  the  i  nf  o  rnati  on  of  Table  I  was 
provided  bv  M^.  Bob  Pose  of  Naval  Ocean  Systems  Center 
(■MOSC).   The   diagrams   included   in   this   appendix  are  sore 


examples   of   the   graphical   output    available 


roni 


tne 


ionospheric  prediction  program  PHOPHTT. 

Figures  46  and  49  are  ray  trace  plots  of  tne  signal  path 
from  WV'V  at  Boulder,  Colorado,  to  San  Antonio  Texas.  Tne 
necessary  inputs  to  tne  program  are  the  date,  time  and 
geophysical  data.  The  date'  and  time  v/ere  selected  to 
correspond  with  the  3WRI  5"oaced  loop  antenna  data.  T;ie 
sunspot  number  and  x-ray  flux  vere  determined  from  published 
geophysical  data.  The  12.7  cm  flux  is  a  number  that  can  te 
determined  from  the  sunspot  number.  PROFHST  also  requires  tne 
transmitter's  location,  power  and  gain  which  in  the  case  of 
A'VV  was  determined  from  published  sources.  The  traces 
represent  tne  signal  at  launch  angles  from  Z  to  52  degrees  in 
5  degree  increments*  the  traces  are  framed  in  a  spatial 
coordinate  s'^'stem,  altitude  versus  range.  The  receiver's 
location  is  denoted  bv  an  asterisk  at  the  correct  distance 
along  the  range  axis.  This  distance  is  the  great  circle  ^vc 
connecting  the  transmitter  and  receive^. 
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Refraction  causes  the  traces  to  tend  tack  to  the  Earth 
unless  the  launch  an^le  is  sufficiently  high  to  permit  the 
ray  to  escape  the  Earth.  This  is  the  case  for  tne  50  degree 
launch  ray  in  fi,Nure  46.  Multincde  occurs  when  ray  traces 
cross  each  other  and  return  to  the  surface  at  approxi'^iately 
the  sane  position.  In  figure  49  the^^e  is  a  Dotential 
rultinode  condition  at  the  receive'-.  The  severity  of 
mult  If  ode  interf  e^^ence  depends  or.  the  strengths  of  tne 
nultihop  traces  after  reflectior  at  the  surface  of  the  Earth 
and  D  layer  absorption. 

Fi5:ures  47  and  50  are  relative  po-/<er  diagrams  in  a 
frequency  vers"'is  24  hour  coordinate  svstep.  The  cur-'yes  are 
relative  rover  contours.  The  tOT)  cuve  represents  the  •'iU?,  and 
its  relative  power  value  is  -3^  do.  The  bottom  contour  is  tne 
LU? ,  and  it  too  has  a  relative  power  value  of  -"6?'  db.  as  one 
works  inward  from  the  outer  two  contours,  each  succeeding 
contour  is  a  +1?  db  higber  than  its  outside  neignbor.  '?ne 
inside  contours  nay  reach  values  of  +10  and  -"-SZ  db .  Tne 
actual  power  received  is  dependent  upon  the  transmitter  power 
and  the  propagation  conditions. 

Plots  of  ionospheric  induced  variance  are  shown  in  figures 
48  and  cl  .  Tne  values  of  variance  are  especially  useful  fcr 
direction  finding  work  in  that  they  are  a  measure  of  error 
induced  by  the  ionosT)here  independent  of  any  RED?  system.  Tne 
values   a  "e  calculate'^,  f"on  empirically  derived  formi;la5.  The 
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fi^st  order  effect  is  frequencv  ieperdent,  and  t^ie  secora 
o^der  effect  is  cased  on  the  Prss  curve  of  variance  as  a 
function  of  range.  Reference  13  is  th<^  MOSC  documentation  on 
the  development  of  the  emDi"i3al  formulas .  Jirures  46  and  51 
show  that  variance  is  typicallv  oetveen  1  and  2  degrees 
squared  with  occassional  sharp  peaks  reaching  ?  degrees 
squa^'ed  for  about  one  honr  duration.  The  najor  peaK  at  13'3^Z 
is  due  to  suni'ise  effects  ( te  rnina  t  c  ;  . 

Additional   information  on  the  capabilities  ^f  ?RO?.-iZT  can 
be  obtained  frorr  ^JOSC. 
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